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Abstract 

We study the Poissonnian ensembles of Markov loops, the related occupation 
field and its renormalized powers in relation with the Gaussian free field. 



1 Introduction 



The purpose of these notes is to explore some simple relations between Markovian path 
and loop measures, spanning trees, determinants, and Markov fields such as the free 
field. They include the notes of the lecture I gave in St Flour in July 2008 with some 
additional material. I choose this opportunity to express my thanks to Jean Picard, to 
the audience and to the readers of the preliminary set of notes, in particular to Jiirgen 
Angst, Cedric Bordenave, Cedric Boutiller, Antoine Dahlqvist, Michel Emery, Jacques 
Franchi, Liza Jones, Adrien Kassel, Rick Kenyon, Thierry Levy, Gregorio Moreno, Bruno 
Shapira, Vincent Vigon, Lorenzo Zambotti and Jean Claude Zambrini. 
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The main emphasis is put on the study of occupation fields defined by Poissonian 
ensembles of Markov loops. These were defined in [15] for planar Brownian motion in re- 
lation with SLE processes and in |16] for simple random walks. They appeared informally 
already in [SH]. For half integral values | of the intensity parameter a, these occupation 
fields can be identified with the sum of squares of k copies of the associated free field (i.e. 
the Gaussian field whose covariance is given by the Green function). This is related to 
Dynkin's isomorphism (cf [7], [26], [H]). We first present the results in the elementary 
framework of symmetric Markov chains on a finite space, proving also several interesting 
results such as the relation between loop ensembles and spanning trees and the refiection 
positivity property. Then we show that some results can be extended to more general 
Markov processes. There are no essential difficulties when points are not polar but other 
cases are more problematic. As for the square of the free field, cases for which the Green 
function is Hilbert Schmidt such as those corresponding to two and three dimensional 
Brownian motion can be dealt with through appropriate renormalization. 

We will show that the renormalised powers of the occupation field (i.e. the self inter- 
section local times of the loop ensemble) converge in the two dimensional case and that 
they can be identified with higher even Wick powers of the free field when a is a half 
integer. 

2 Symmetric Markov processes on finite spaces 

Notations: Functions and measures on finite (or countable) spaces are often denoted as 
vectors and covectors. 

The multiplication operator defined by a function / acting on functions or on measures 
is in general simply denoted by /, but sometimes it will be denoted Mf. The function 
obtained as the density of a measure fj, with respect to some other measure u is simply 
denoted -. 

2.1 Graphs 

Our basic object will be a finite space X and a set of non negative conductances Cx,y = 
Cy^xi indexed by pairs of distinct points of X. This situation allows to define a kind of 
discrete topology and geometry we will briefiy study in this section and in the following 
ones. 
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We say {x, y} is a link or an edge iff C^^y > and an oriented edge (x, y) is defined by 
the choice of an ordering in an edge. We set —{x,y) = {y,x) and if e = {x,y), we denote 
it also (e~, e"*"). The degree of a vertex x is by definition the number of edges incident 
at X. 

The points of X together with the set of non oriented edges E define a graph (X, E) . 
We assume it is connected. The set of oriented edges is denoted E°. It will always be 
viewed as a subset of X^, without reference to any imbedding. 

The associated line graph is the oriented graph defined by E° as set of vertices and 
in which oriented edges are pairs (ei, 62) such that = 63 • The mapping e — — e is an 
involution of the line graph. 

An important example is the case in which conductances are equal to zero or one. 
Then the conductance matrix is the adjacency matrix of the graph: C^^y = l{x,y}eE 

A complete graph is defined by all conductances equal to one. 

The complete graph with n vertices is denoted i^„.The complete graph K4 is the graph 
defined by the tetrahedron. is not planar (i.e. cannot be imbedded in a plane), but 
is. 




A finite discrete path on X, say {xo,Xi, ...Xn) is called a (discrete) geodesic arc iff 
{xi,Xi+i} G E (path segment on the graph) and 7^ Xj+i (without backtraking) . 
Geodesic arcs starting at Xq form a marked tree %xq rooted in xq (The marks are the 
points of X). Oriented edges are defined by pairs of geodesic arcs of the form: 

((xo, xi, ...x„), (xo, xi, ...x„, x„+i)) (the orientation is defined in reference to the root). 
T^n is a universal cover of X [25l. 
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On the space of discrete loops based at some point Xq, we can define an operation 
of concatenation, which provides a monoid structure, i.e. is associative with a neutral 
element (the empty loop). The concatenation of two closed geodesies based at xq is 
not directly a closed geodesic. It can involve backtracking "in the middle" but then after 
cancellation of the two inverse subarcs, we get a closed geodesic, possibly empty if the two 
closed geodesies are identical up to reverse order. With this operation, closed geodesies 
based at Xq define a group Fa-^.The structure of does not depend on the base point and 
defines the fundamental group T (as the graph is connected: see for exemple [25]). Any 
geodesic arc 71 from xq to another point of X defines an isomorphism between T^q and 
Ty^. It associates to a closed geodesic 7 based in xq the closed geodesic [7i]~^77i (here 
[71]"^ denotes the backward arc. In the case where xq = yo, it is an interior isomorphism 
(conjugation by 71). 

There is a natural left action of T^^ on T^^g . It can be interpreted as a change of root 
in the tree (with the same mark). Besides, any geodesic arc between xq and another 
point yo of X defines an isomorphism between T^q and Tyg (change of root, with possibly 
different marks) . 

We have just seen that the universal covering of the finite graph {X, E) at Xq is a tree 
%xo projecting on X. The fiber at xq is Yxq- These groups are conjugated in a non canonical 
way. Note that X = Tj-^yi^g (here the use of the quotient on the left corresponds to the 
left action). 

Example 1 Among graphs, the simplest ones are r— regular graphs, in which each point 
has r neighbours. A universal covering of any r— regular graph is isomorphic to the 
r— regular tree X'-'"-'. 

Example 2 Cayley graphs: A finite group with a set of generators S = {gi,..gk} such 
that S n is empty yields an oriented 2k-regular graph . 

A spanning tree T is by definition a subgraph of {X, E) which is a tree and covers all 
points in X. It has necessarily |X| — 1 edges. See for example two spanning trees of K4. 
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Two spanning trees of 



Its inverse images by the canonical projection from a universal cover Tx„ onto X form 
a tesselation on X^^, i.e. a partition of in identical subtrees, which are fundamental 
domains for the action of T^q- Conversely, a section of the canonical projection from the 
universal cover defines a spanning tree. 

Fixing a spanning tree determines a unique geodesic between two points of X. There- 
fore, it determines the conjugation isomorphisms between the various groups r^-g and the 
isomorphisms between the universal covers %xo- 

The fundamental group T is a free group with — |X| + 1 = r generators. To 
construct a set of generators, one considers a spanning tree T of the graph, and choose an 
orientation on each of the r remaining links. This defines r oriented cycles on the graph 
and a system of r generators for the fundamental group. (See [25] or Serres ([36]) in a 
more general context). 

Example 3 Consider and K4^. 

Here is a picture of the universal covering of K4, and of the action of the fundamental 
group with the tesselation defined by a spanning tree. 
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Universal cover and tesselation of K4 



There are various non ramified coverings, intermediate beween {X, E) and the uni- 
versal covering. Non ramified means that locally, the covering space is identical to the 
graph (same incident edges). Then each oriented path segment on X can be lifted to the 
covering in a unique way, given a lift of its starting point. 
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Each non ramified covering is (up to an isomorphism) associated with a subgroup H of 
r, defined up to conjugation. More precisely, given a non ramified covering X, a point Xq 
of X and a point xq in the fiber above Xq, the closed geodesies based at Xq whose lift to the 
covering starting at Xq are closed form a subgroup ifjj of F^^, canonicaly isomorphic to 
the fundamental group of X represented by closed geodesies based xq. If we consider a 
different point ijQ, any geodesic path segment 71 between xq and ^ defines an isomorphism 
between r^.^ and Vy^ which exchanges and Hy^. Denoting 71 the projection of 71 on 
X , it associates to a closed geodesic 7 based in xq whose hft to the covering is closed the 
closed geodesic [71] "^771 whose hft to the covering is also closed. 

Conversely, if if is a subgroup of Vx^i the covering is defined as the quotient graph 
{Y, F) with Y = if\Tx„ and F the set of edges defined by the canonical projection from 
oiito Y. H can be interpreted as the group of closed geodesies on the quotient graph, 
based at -f^^oj i-e. as the fundamental group of Y. 

If if is a normal subgroup, the quotient group (also called the covering group) ifyFj-g 
acts faithfully on the fiber at xq. An exemple is the commutator subgroup [Txq, '^xq]- The 
associate covering is the maximal Abelian covering at xq 

Example 4 By central symmetry, the cube is a two fold covering of the tetrahedron as- 
sociated with the group Z/2Z. 

Exercise 5 Determine the maximal Abelian cover of the tetrahedron. 
2.2 Energy 

Let us consider a nonnegative function k on X. Set = + C^^y and Py = P 
is a A-symmetric (sub) stochastic transition matrix: XxPy = with P^ = for all x 

in X. 

It defines a symmetric irreducible Markov chain 

We can define above it a continuous time A-symmetric irreducible Markov chain Xf, 
with exponential holding times of parameter 1. We have Xf — ^Nt, where denotes a 
Poisson process of intensity l.The infinitesimal generator is given by = Py — 6y. 

We denote by Pf its (sub) Markovian semigroup exp{Lt) — |t-^ ■ L and Pt are 
A-symmetric. 

We will use the Markov chain associated with C, k, sometimes in discrete time, some- 
times in continuous time (with exponential holding times). 
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Recall that for any complex function G X, the "energy" 
is nonnegative as it can be written (easy exercise) 

x,y X X x,y 

The Dirichlet space ([S]) is the space of real functions equipped with the energy scalar 
product 

x,y X x x,y 

defined by polarization of e. 

Note that the non negative symmetric "conductance matrix" C and the non negative 
equilibrium or "killing" (or "equilibrium" ) measure n are the free parameters of the model. 
The eigenfunction associated with the lowest eigenvalue of —L has constant sign by the 
well known argument based on the fact that the map z z'^ lowers the energy. 

We have a dichotomy between: 

- the recurrent case where is the lowest eigenvalue of — L, and the corresponding 
eigenspace is formed by constants. Equivalently, PI = 1 and k, vanishes. 

- the transient case where the lowest eigenvalue is positive which means there is a 
"Poincare inequality": For some positive e, the energy e{f,f) dominates e{f,f)^ 
for all /. Equivalently, (as we are on a finite space) k does not vanish. 

In the transient case, we denote by V the associated potential operator (— L)^^ = 
Jq°° Ptdt. It can be expressed in terms of the spectral resolution of L. 

We denote by G the Green function defined on as = ^ = ^[(/ — Py^]y 
i.e. G = (Ma — C)^^. It induces a linear bijection from measures into functions. We will 
denote "^yG^'^^y by {GuY or G^{x). 

Note that Gfi is characterized by the identity e{f,Gfi) = (/, /i) (i.e. J2xf^f^x) valid 
for all functions / and measures fi. In particular Gn = 1 as e(l, /) = X] f^'^x = (/? 1)^- 
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See ([H]) for a development of this theory in a more general setting. 

In the recurrent case, the potential operator V operates on the space A"*- of functions / 
such that (/, 1)^ = as the inverse of the restriction of J — P to A"*". The Green operator G 
maps the space of measures of total charge zero onto A"*": Setting for any signed measure 
u of total charge zero Gu = Vj, we have for any function /, {u, f) = eiGv, f) (as 
e(Gz/, 1) = 0) and in particular — P = e{G{5x — Sy), /). 

Exercise 6 Compute the Green operator in the case of the complete graph 

In quantum mechanics, —L is called the Hamiltonian and its eigenvalues are the energy 
levels. 

One can learn more on graphs and eigenvalues in [2] 
2.3 Feynman-Kac formula 

A discrete analogue of the Feynman-Kac formula can be given as follows: Let s be any 
function on X taking values in (0, 1]. Then, for the discrete Markov chain C,n associated 
with P, it is a straightforward consequence of the Markov property that: 

n-l 

E,{l[s{Qi{^^=y}) = [{M,prry 

j=0 

Similarly, for the continuous time Markov chain xt (with exponential holding times), if 
k{x) is a nonegative function defined on X, we have: 

E^{e-foH^^)<in{^^^yy) = [exp(t(L - M^)]^. 

It can be shown easily by differentiating the first member V{t) with respect to t, to check 

V'{t) = {L-Mk)V{t). 

For any nonnegative measure set = (— L + and G^ = V^Mi = {M\ + 

— G)^^. It is a symmetric nonnegative function on X x X. Gq is the Green function 

G, and can be viewed as the Green function of the energy form = e + || IIl2(^). 
Note that has the same conductances G as e, but x is added to the killing measure. 

Note also that is not the potential of the Markov chain associated with when one 

takes exponential holding times of parameter 1: the holding time parameter at x becomes 
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1 + x(^)- But the Green function is intrinsic i.e. invariant under a change of time scale. 
Still, we have by Feynman Kac formula 

/"OO 

^0 

We have also the "resolvent" equation V -V^^ VMkV^ = V^MkV. Then, 

G-G^ = GM^G^ = G^M^G (1) 

Note that the recurrent Green operator G defined on signed measures of zero charge 
is the limit of the transient Green operator G^, as x — > 0. 

2.4 Recurrent extension of a transient chain 

It will be convenient to add a cemetery point A to X, and extend C, A and G to — 
{X U A} by setting , A a = Exex '^^^ = and G^'^ = G^'^ = G^'^ = for all 
xeX. Note that \{X^) = Exxx Cx,y + 2 Ex i^x) 

One can consider the recurrent "resurrected" Markov chain defined by the extensions 
the conductances to X^. An energy is defined by the formula 



2 

From the irreducibility assumption, if follows that vanishes only on constants. We 
denote by the transition kernel on X^ defined by 

or equivalently by 

rpAla; _ ^x,y _ Cx,y 

Z^yeX^ ^x,y 

Note that P^l = 1 so that A is now an invariant measure.with Xx[P^]y = \[P^]l on X^ 
Denote and G^ the associated potential and Green operators. 

Note that for /i carried by X, for all x e X, denoting by £a the unit point mass at A, 

/X. = e^(G'^(^ - /i(X)£A), 1.) = A,((/ - P^)G'^(/x - fi{X)eA)ix) 
= Xx{{I - P)G'^{fi - fi{X)eA)){x) - KxG^ti - fi{X)eA){A). 
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Hence, applying G , it follows that on X^, 

Gfi = G^{n-n{X)6A)-G^il^-f^iX)eA){A)GK = G^{fi-fi{X)6A)-G^{fi-n{X)eA)iA). 
Moreover, as G^{fi — fi{X)eA) is in A"*-, 

J2 AxG(/i)" = -G^ifi - /i(X)£A)(A)A(X^). 

Therefore, G'^ifi - /i(X)£A)(A) = ^gg^ and - fi{X)eA) = + Gfi. 

This type of extension can be done in a more general context ( See [21j and Dellacherie- 
Meyer [5]) 

Remark 7 Conversely, a recurrent chain can be killed at any point xq of X , defining a 
Green function G^~^^°^ on X — {xq}. Then, for any /i carried by X — {xq}, 

G^-{-o}^ = ^(/i - - Gifi - /i(X)£,.J(xo). 

This transient chain allows to recover the recurrent one by the above procedure. 

Exercise 8 Consider a transient process which is killed with probability p at each passage 
in A. Determine the associated energy and Green operator. 

2.5 Transfer matrix 

Let us suppose we are in the recurrent case: We can define a scalar product on the space 
A of functions on E° (oriented edges) as follows 

{uj,!])^ = lJ2xy^'->:,y^^'^''l^'^- Denoting as in [23] df^''" = f" — we note that 
{df,dg)^ = e{f,g) In particular 

{df,dGi5y-s,))^ = df^^y 

Denote A_, (A_|_) the space of real valued functions on E° odd (even) for orientation 
reversal. Note that the spaces and A_ are orthogonal for the scalar product defined 
on A. The space A_ should be viewed as the space of "discrete differential forms". 
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For any a. in A_, define d*OL by (d*Q;)^ = ~^y^x^y^^''^ ■ Note it belongs to A"*" as 
Yl,x,y Cx,ya'''^ vanishes. 
We have 

xeX x,y xex 

as the two terms of the sum are in fact equal since a is skew symmetric. The image of 
d and the kernel of d*are therefore orthogonal in A_. We say a in A_ is harmonic iff 
d*a = 0. 

Moreover, e(/, /) = {df,df)^ = E.ex(^*^/)V"A,. 
Note that for any function f , d*df ^ -Pf + f ^ —Lf . 

The projection of any a in A_ on the image of d, is easily obtained as dVd*{a). Indeed, 
for any function g, {a^dg)^ = {d*a,g)-^ = e{Vd*a,g) = {dVd*{a),dg)j^. 
d is the discrete analogue of the differential and d* the analogue of its adjoint, depending 
on the metric which is here defined by the conductances. 

Set CK^u^^) = {^tV) — ='=('^)t^) ^-nd elsewhere. Then Xj;d*a(^u^y){x) — — 6^ 

and dVd*{a(^u,v)) = dG{S^ — Note that given any orientation of the graph, the family 
{'^(uv) ~ \/^^(^iu,v), (u.v) G E^} is an orthonormal basis of A_ (here denotes the 
set of positively oriented edges). 

The symmetric transfer matrix K^^'y'>'^^''"\ indexed by pairs of oriented edges, is defined 
to be 

j^{x,yUu,v) ^ [^^(^^ _ ^^-^py ^ ^(^^ _ ^^y _ q^^^ _ ^^(^^ _ ^^(^^ _ 

for x,y,u,v G X, with Cx,yCu,v > 0. For every oriented edge e = (a;, in X, set 
= dG{5y -d""). 

We have {K^i K^)^ — K^'^.We can view dG as a linear operator mapping the space 
measures of total charge zero into A_. As measures of the form 5y — 5x span the space of 
measures of total charge zero, it is determined by the transfer matrix. 

Note that d*dGv — v/X for any v of total charge zero, that for all a in A_, {d*a)X 
has total charge zero and that dG{{d*a)X) — dVd*{a) is the projection Il{a) of a on the 
image of d in A_. 
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In particular 



< n(a(„,y)) >A=< a^u,v), dG{6y> - 6u') >a= 

Consider now, in the transient case, the transfer matrix associated with C^. 
We see that for x and y in X, G^{6^ - 5^,)" - - dyY = G(4 - 5^,)" - G(4 - Sy)". 

We can see also that G"^(5x — 5a) = G6x — ^^(xA)'' • So the same identity holds in X'^. 
Therefore, as G^'^ = 0, in all cases, 

j^ix,y),{u,v) _ QX,u _|_ Qy,v _ qx,v _ Qy,u 

Exercise 9 Cohomology and complex transition matrices. 

Consider, in the recurrent case, u G A_ such that d*uj = 0. Note that the space of 
such uj's is isomorphic to the first cohomology space, i.e. the quotient A_/Im((i). Prove 
that P{I + iuj) is X—self adjoint on X , maps 1 onto 1 and that we have ^xiYYjZo 1 + 
u;(^j,^j+i)liU=y})^mi + iu;))X 

3 Loop measures 

3.1 A measure on based loops 

We denote by the family of probabihty laws on piecewise constant paths defined by 
Pf 

^""(7(^1) = ^1, ■■■,l{th) = Xh) = Pt^{x,Xi)Pt^-ti{xi,X2) . . . Ptf^-tf^_,{Xh-l, Xh) 

Denoting by ^(7) the number of jumps and Tj the jump times, we have: 
^xivil) = k, -fTi = xi, -fTk-i = Xk-1, Ti e dti, Tk e dtk) 

■^{o<h<...<tk}^~*''dti...dtk 



^X^Xl ■■■^Xk 



For any integer p > 2, let us define a based loop with p points in X as a couple I — {^,t) — 
((Cm, 1 < ?7i < p), (Tm, 1 < m < p + 1), ) in XP X M.^^, and set = ^p+i (equivalently, we 
can parametrize the associated discrete based loop by Z/pZ). The integer p represents 
the number of points in the discrete based loop ^ = (d, ■■■^p{^)) and will be denoted p{^). 
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Note two time parameters are attached to the base point since the based loops do not in 
general end or start with a jump. 

Based loops with one point (p = 1) are simply given by a pair {C,, t) in X x M_|_. 

Based loops have a natural time parametrization l{t) and a time period T(^) = 
^fif^Vi. If we denote YliLin by l{t) = ^m-i on [T^_i,r^) (with by convention 
To = and = Q- 

Let P^'^ denote the (non normahzed) "bridge measure" on piecewise constant paths 
from X to y of duration t constructed as follows: 

If ti < t2... <th<t, 

Ay 

Its mass is p^'^ = ~\- measurable set A of piecewise constant paths indexed by 

[0 t] , we can also write 

Ay 

A cr-finite measure /i is defined on based loops by 

Prom the expression of the bridge measure, we see that by definition of /i, if ti < 
t2... <th<t, 

= x^, ^x,,Te dt) = [Pi,+*_*J^jP*,_iJ^^...[Pi,_*,_J^;;-^di. (2) 

Note also that for A; > 1, using the second expression of P^'^ and the fact that condition- 
ally on Nf = k, the jump times are distributed like an increasingly reordered /c— uniform 
sample of [0 t] 

A^P,^'^(p = A;, 6 = X2, 4 = Xk, Ti e dh, Tk e dh) 

^e-'-p:p:^...p:n^o<H<..,,<t}j^dt,...dt, 
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Therefore, 



= k,^i = Xi, ..,^k = Xk,Ti e dti, .., Tfc G dtk, T G dt) 



(3) 




(4) 



for k > 1. 

Moreover, for one point-loops, = Ij^i = ^i, £ dt} = ^dt. 



It is clear on these formulas that for any positive constant c, the energy forms e and 
ce define the same loop measure. 

3.2 First properties 

Note that the loop measure is invariant under time reversal. 

If D is a subset of X, the restriction of fi to loops contained in D, denoted /i^ is 
clearly the loop measure induced by the Markov chain killed at the exit of D. This can 
be called the restriction property. 

Let us recall that this killed Markov chain is defined by the restriction of A to D and 
the restriction of P to (or equivalently by the restriction cc of the Dirichlet norm 
e to functions vanishing outside D). 

As J ^^e~^dt = i, it follows from ([3]) that for /c > 1, on based loops. 



Mi) = k,ii = = Xk) = ^Pxl--P, 



->Xk 



(5) 



In particular, we obtain that, for A; > 2 



^{p = k) = \Tr{P'') 



and therefore, as Tr{P) = 0, 



fxip >l) = Yl vTriP') = -log(det(/ - P)) = log(det(G) J] A.) 



(6) 



2 X 



since (denoting Mx the diagonal matrix with entries A^.), we have 



det(/-P) 



det(MA - C) 
det{Mx) 



15 



Note that det{G) is defined as the determinant of the matrix G^'"^. It is the determinant 
of the matrix representing the scalar product defined on M'"^' (more precisely, on the space 
of measures on X) by G in any basis, orthonormal with respect to the natural euclidean 
scalar product on M'"^'. 

Moreover 

/oo 
p(/)l|,>i}/x(rfO = J^TriP') = Tr((/ - P)-'P) = Tr{GC) 
2 

3.3 Loops and pointed loops 

It is clear on formula [2] that /i is invariant under the time shift that acts naturally on 
based loops. 

A loop is defined as an equivalence class of based loops for this shift. Therefore, /x 
induces a measure on loops also denoted by /i. 

A loop is defined by the discrete loop ^° formed by the in circular order, (i.e. up to 
translation) and the associated holding times. We clearly have: 



Mr = (a:i,x2,...,x,)°) = p;^..p;^'= 



However, loops are not easy to parametrize, that is why we will work mostly with 
based loops or pointed loops. These are defined as based loops ending with a jump, or as 
loops with a starting point. They can be parametrized by a based discrete loop and by 
the holding times at each point. Calculations are easier if we work with based or pointed 
loops, even though we will deal only with functions independent of the base point. 

The parameters of the pointed loop naturally associated with a based loop are ^i, 
and 

n + rp+i= T*,Ti = T*, 2<i<p 

An elementary change of variables, shows the expression of /i on pointed loops can be 
written: 

^{p = k,^, = x„r; G dt,) = P:i...P:-^^e-^'^dt,...dtk. (7) 
Trivial (p = 1) pointed loops and trivial based loops coincide. 
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Note that loop functionals can be written 

$(/°) = ^l{^=,}$,((e„T;),i-l,...A;) 

with invariant under circular permutation of the variables (^j,Tj*). 
Then, for non negative 

and by invariance under circular permutation, the term ti can be replaced by any ti. 
Therefore, adding up and dividing by k, we get that 

J Mf)ii{di) = j^-Y,^k{xi,u)p:^...p:^e-^'^dti...dh. 

The expression on the right side, applied to any pointed loop functional defines a 
different measure on pointed loops, we will denote by It induces the same measure as 
// on loops. 

We see on this expression that conditionally on the discrete loop, the holding times of 
the loop are independent exponential variables. 

/(p = Ki, = Xi, T* e dt,) = \\{ ^^^e-'^dU (8) 



Conditionally on p(^) = fc, T is a gamma variable of density ^*_j^^, e * on IR+ and 

( j^, 1 < i < /c) an independent ordered /c-sample of the uniform distribution on (0,T) 
(whence the factor |). Both are independent, conditionally on p of the discrete loop. We 
see that /i, on based loops, is obtained from /i on the loops by choosing the based point 
uniformly. On the other hand, it induces a choice of biased by the size of the t*'s, 
different of pL* (whence the factor |). But we will consider only loop functionals for which 
H and yU* coincide. 

It will be convenient to rescale the holding time at each by A^. and set 
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The discrete part of the loop is the most important, though we will see that to estab- 
lish a connection with Gaussian fields it is necessary to consider occupation times. The 
simplest variables are the number of jumps from x to y, defined for every oriented edge 

Nx,y = #{i ■ii = x, ii+i = y} 
(recall the convention ^p+i = ^i) and 

y 

Note that = > ^ '■ = x} except for trivial one point loops for which it vanishes. 
Then, the measure on pointed loops ([7]) can be rewritten as: 

dt 

= 1,^ = x,f e c/t) = e"^"*y and (9) 

Another bridge measure fi^'^ can be defined on paths 7 from x to y: fi^'^ld'y) = 
J,^n''{d^)dt. 

Note that the mass of /x^'^ is G^'"^. We also have, with similar notations as the one defined 
for loops, p denoting the number of jumps 

/i^'f (p(^) = k, = Xi, 7Tfc_i = Xk-i, Ti G rfti, Tk-i E dtk-i,T G dt) 

Cx,X2Cx2,X3---Cxf,_^,y _i 

— r-^ r ^{o<ti<...<tf,<t}G ati...atkat. 

'^x'^X2 • • • '^y 

For any x ^ y m X and s G [0, 1], setting pi = P" if (m, v) ^ {x, y) and Py^^'^ = sPy, 
we can prove in the same way as that: 

Ms'^-n|,>i}) = -log(det(/-P(^))). 

Differentiating in s = 1, and remenbering that for any invertible matrix function M{s), 
A 

ds 



log(det(M(s)) = Tr{M'{s)M{s)-^), it follows that: 



^i{Nx,y) = [(/ - P)-']lPy^ = G^^yCx,y 
and fi{Nx) = Eyf^iNx,y) = A.G-'^ - 1 (as G{Mx - C) = Id). 
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Exercise 10 Show that more generally iJt,{N^^y{N^^y—l)...{N^^y—k-\-l)) = {k—l)\{G^''^Cx,yY 
Exercise 11 Show that more generally, if Xi,yi are n distinct oriented edges: 



We finally note that if Cx,y > 0, any path segment on the graph starting at x and 
ending at y can be naturally extended into a loop by adding a jump from y to x. We have 
the following 

Proposition 12 a) For Cx,y > 0, the natural extension of /^^'^ to loops coincides with 
^^^aidl). 

b)More generally, ifxi,yi are n oriented edges, Y\ '^Z^'^'^''^ fi{dl) can be obtained as the 
sum of the images, by concatenation in all circular orders, of the product of the bridge 
measures 

Proof. The first assertion follows from the formulas, noticing that a loop / can be 
associated to Ny^^iP) distinct bridges from x to obtained by "cutting" one jump from 
y to X. ■ 

3.4 Occupation field 

To each loop f we associate local times, i.e. an occupation field {IxiX e X} defined by 

-.T(0 -, p(0 



Jo his) ^ 



for any representative / = (^j,T*) of 1°. 

For a path 7, 7 is defined in the same way. 
Note that 

M(l-e-"'^)l{,=i})= re-*(l-e-^*)y = log(l + f ) (H) 

Jo -^x 

(by expanding 1 — e before the integration, assuming first that a is small and then 
by analyticity of both members, or more elegantly, noticing that J^(e''^ — e"''^)^ is 
symmetric in (a.b) and (c.d)). 
In particular, /x(Pl{p=i}) = j-. 
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From formula [3, we get easily that the joint conditional distribution of (/^, x G X) 
given [Nx, x G X) is a product of gamma distributions. In particular, from the expression 
of the moments of a gamma distribution, we get that for any function $ of the discrete 
loop and k > 1, 

/i((r)'=l|p>i}$) = X~'fi{{N, + k- 1)...{N, + l)iV,$) 

In particular, = -^[fiiN^) + 1] = G^'"^. 

Note that functions of I are not the only functions naturally defined on the loops. 
Other such variables of interest are, for n >2, the multiple local times, defined as follows: 

- ""W 1 

= ^ / i{ati)=xr+^,....at„.j)=x„,...atr,)=x,} JJ -^^u 

j=0 J0<ti<...<t„<T ^x. 

It is easy to check that, when the points Xi are distinct, 

n— 1 n 

^"'■■■"" = E E UM,n-^...}nr (12) 

j=0 l<ii<..<i„<p{l) 1=1 

Note that in general /^i'-- '^fc cannot be expressed in terms of /. 

If Xi = X2 = . . . = Xn, /^i'---'^" = j-^^j^^i [/^]". It can be viewed as a n-th self intersection 
local time. 

One can deduce from the defintions of fi the following: 

Proposition 13 = G'^i.^'2G'^2,x3_ 

In particular, = i ^^^^^ G'^-(i)'^-(2)G^-(2)'^-(3)...G^^(n),^.(i) 

Proof. Let us denote ^ [Pt\y by p^'^ or pt{x,y). From the definition of /^i'----^" and 
/X, equals: 

n-1 

/ / 



a: j=Q 



{0<ti. ..<*„<*} 



where sums of indices k + j are computed mod(n). By the semigroup property, it equals 

X] / / ^Pt2-tAxi+j,X2+j) . . .pt,+t-t„{Xn+j,Xi+j)Y[dUdt. 



j=0 



{0<ti<...<t„<t} t 
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Performing the change of variables V2 = t2 — ti, ..,Vn = tn — t"! = ti + t — t„, and 
V = ti, we obtain: 

+i) I I dvidv 

~^J{0<v<vi,0<v,}'Vl + ---+Vn 
n—1 „ 

= \ ~ P^2 iXl+j,X2+j) Pvi (Xn+j, Xi+j) TT dVi 

~^J{0<v,}Vl + ■■.+Vn 

n „ 

= / ■ ^— p^^{Xi,X2) . . .PvAXn,Xl)Y\d'Vi 

= Pv2{Xl,X2) . . .Pvi{Xn,Xi)Y\dVi 

J{0<Vi} 

Note that another proof can be derived from formula (1121) ■ 

Exercise 14 (Shuffle product) Given two positive integers n > k, let Vn,k be the family 
of partitions of into k consecutive non empty intervals Ii = {ii,ii+i, ...,ii^i — 1) with 
ii = 1 < i2... <ik < ik+i =n + l. 
Show that 

m-l inf(n,m) 

'jxi,...,xn'jyi,---,ym _ 'l^ii'yj+Ji^^i2'-yi+Jk 
j=o k=i ier„^k Je-P,n,fc 

^n-l 



Similarly, we can define A^(^i,y,),...(x„,y„) to be Y.j=o Ei<n<..<i„<p(0 Uti ^{i.,=x,+„i.^+^=m+j}- 
If = {x,y) for all ^, it equals ^-.^(^-.^-i^-^jf^-.-^^+i) 

Notice that _^ 



Then we have the following: 
Proposition 15 / iV(.„,,),...(.„,,„, = (11 a.^JG^^'^^G^^'^^-G^" 

The proof is left as exercise. 
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Exercise 16 For Xi — X2 — ■■■ — x^, we could define different self intersection local times 

k 



l<ii<..<ijfc<p(Z) 1=1 



which vanish on < k. Note that 



i=l 



1. For any function $ of the discrete loop, show that 

2. More generally prove in a similar way that 

Hi^L '(p) _ A3. n{ — i{jv,>fc}5'j- 

Let us come back to the occupation field to compute its Laplace transform. Prom 
the Feynman-Kac formula, it comes easily that, denoting Mx the diagonal matrix with 
coefficients ^ 

Pr(e-(^>^) - 1) = f (exp(i(P - / - MJ)^ - exp(i(P - /))^). 

Ax ^ 

Integrating in t after expanding, we get from the definition of /j, (first for x small enough): 

(e-(''>^) - l)dMO = E / [Tr{{P - M^r) -Tr{{Pn'-—e~'dt 
fc=i -^0 ^ 

CO ^ 

fe=i 

= -Tr(log(7 - P + M^)) + rr(log(7 - P)) 

Hence, as T'r(log) = log(det) 

J (e"('^^) - l)dii{l) = log[det(-L(-L + M^/a)"^)] = -logdet(/ + VM^) 
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which now holds for all non negative x ^'^ both members are analytic in x- Besides, by 
the "resolvent" equation ([1]): 



det(/ + GMj-i = det(/ - G,M,) = (13) 

Note that det(/+GM^) = det (J+M^^GM^^) and det(/-G^M;^) = det(/-M^G'^My^), 
so we can deal with symmetric matrices. Finally we have 



Proposition 17 /i(e~(''^) - 1) = -log(det(J + M^GM^)) = log( 



det(Gx)' 



det(G) ' 

Note that in particular /i(e~*''" — 1) = — log(l + tG^'^). Consequently, the image 



det(Gf^ 



measure of n by is l{s>o}^-^—ds. 

Note finally that if x has support in D, by the restriction property 

/^(l{r(x\i.)=o}(e"^'"^ - 1)) = -log(det(/ + M^G^M^)) = log(^^^). 

Here the determinants are taken on matrices indexed by D and G'^ denotes the Green 
function of the process killed on leaving D. 



For paths we have P^''^(e"V''^)) = ^ exp(t(L - M^))^^y. Hence 

^x,.(,-(7,x>) = 1 ((J _ p + M^/,r%,y = [G^py. 

Ay 

In particular, note that from the resolvent equation ([1]), we get that 

Hence ^^yj^ = ^^^Qa:,x and under the probabihty ^5^, Ix follows an exponential 
distribution of mean G^'^. 

Also E^(e-<^'^>) = 'EylGxT'^f^y i.e.[G^K]^. 

Finally, let us note that a direct calculation shows the following 

Proposition 18 On loops based in x, j^^'^{dl) = l^iJ,{dl) 

More generally, if Xi are n points, F^''"''^"fi{dl) can be obtained as the the image by 
circular concatenation of the product of the bridge measures fi^^'^^+^{dl) and J^P'//((iZ) 
can be obtained as the sum of the images, by concatenation in all circular orders, of the 
product of the excursions measures 
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Remark 19 Propositions and can be generalized: For example, one can evaluate 
expressions of the form Y[ 11 -rr^^^ l^{dl) as a sum of images, by concatenation in all 
circular orders, of a product of bridge measures . 



3.5 Wreath products 

Associate to each point of X an integer n^. Let Z be the product of the groups Z/riajZ. 
On X X Z, define a set of conductances C(^x,z),{x',z') by: 



C(x,z),{x',z') - ——^x,x' J_J_ ^{Zy=Z^} 

' I'X ' I'x' , , 

y^x,x' 

and set li[x,z) = i^x- Let e be the corresponding energy form. Note that \x,z) = ^x- 
Then, denoting Ji the loop measure and P the transition matrix on X x Z defined by 
e, we have the following 

Proposition 20 / l{p>i} nx,7V4z)>o ^l^i^^) = J^(P > 1) = -log(det(/ - P)). In partic- 
ular, if = n for all x, 

J l{p>i}n-#^^' ^-«>oV(rf/) =Jl{p>l) = - log(det(/ - P)) 

Proof. Each time the Markov chain on X x Z defined by e jumps from a point above 
X to a point above y, and Zy are resampled according to the uniform distribution on 

Z/n.Z X Z/uyZ. It follows that ET=AP%i = Px.Px^ -P^'-' n,e{x,x„....,.,} ^- 

The detail of the proof is left as an exercise. ■ 

This construction therefore gives an interesting information about the number of dis- 
tinct points visited by the loop, which is more difficult to evaluate than the occupation 
measure. 

In the case where X is a group and P defines a random walk, P is associated with a 
random walk on X x Z equipped with its wreath product structure (Cf (30j). 



3.6 Countable spaces 

The assumption of finiteness of X can of course be relaxed. On countable spaces, the 
previous results extend easily under spectral gap conditions. In the transient case, the 
Dirichlet space EI is the space of all functions / with finite energy e(/) which are limits 
in energy norm of functions with finite support. The energy of a measure is defined as 



24 



supjgjj^^^. Finitely supported measures have finite energy. The potential Gjj, is well 
defined in for all finite energy measures /x, by the identity e(/, Gyu) = (/, /i), valid for all 
/ in the Dirichlet space. The energy of the measure equals e{Gfi) = {Gfi,fi). 

Most important examples are the non ramified covering of finite graphs (Recall that 
non ramified means that the projection is locally one to one, i.e. that the projection on 
X of each vertex v of the covering space has the same number of incident edges as v 
) . Consider a non ramified covering graph {Y, F) defined by a normal subgroup H^q of 
r^..^,. The conductances C and the measure A can be lifted in an obvious way to Y as 
-f^xo\rxo-periodic functions but the associated Green function G or semigroup are non 
trivial. 

We have G^'^ = ^^gj:/^ (7«(^)'7(«(y)) for any section i of the canonical projection 
from Y onto X. 

Let us consider the universal covering (then if^.^ is trivial). It is easy to check it will 
be transient even in the recurrent case as soon as {X, E) is not circular. 

The expression of the Green function G on a universal covering can be given exactly 
when it is a regular tree, i.e. in the regular graph case. In fact a more general result can 
be proved as follows: 

Given a graph (X, i?), set dx = J2y^{x,y}eE (degree or valency of the vertex x), 
Dx,y = dx5x,y and denote Ax^y the incidence matrix lE{{x,y}). 

Consider the Green function associated with A^^ = {dx — 1)u + ^^ with < m < inf(^^-^) 
and for {x,y} e E, Cx,y = 1: Then, noting that J^x ^z,xu'^'~'''^^ = (4-l)n'^(^'^)+^+M'^(^'^)-\ 
we have G = [u~^I + {D- I)u - A]-\ 

Proposition 21 On the universal covering %xo, G ^''^ = u'^(^'S^)(i — u)^^ . 

Indeed, ^^(Az^^ — Az^x)u'^^^''^^ = for z y and equals ^ — u for z = y. 

It follows that for any section i of the canonical projection from %xo onto X, 

3.7 Zeta functions for discrete loops 

We present briefiy the terminology of symbolic dynamics (see for example pHj) in this 
simple framework (setting /(xq, Xi..., = \og{PxQ^xi)i P induces the Ruelle operator 
Lf associated with /). 
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The pressure is defined as the logarithm of the highest eigenvalue (3 oi P. It is asso- 
ciated with a unique normalized positive cigenfunction h, by Perron Probenius theorem. 
Note that Ph — (3h implies XhP — (3Xh by duality and that in the recurrent case, the 
pressure vanishes and h = —7^ — . 

In continuous time, the lowest eigenvalue of —L i.e. l—f3 plays the role of the pressure 
The equilibrium measure associated with /, m = h'^X is the law of the stationnary 

Markov chain associated with transition probability -^Pyhy. 

If PI = 1, i.e. K = 0, we can consider a Feynman Kac type perturbation P^^*^) = 

PM A , with £ I and k a positive measure. Perturbation theory shows that (5^^'^^ — 1 — 

j^eI - 1 + o{e) = + o{e) and that = ^ + o{e). 

Wc deduce from that the asymptotic behaviour of J {e~'^^'^'^ — = log(dct(/ — 

p(£'^))) _ log(det(/ - p(=(«+x))) which is equivalent to - log(l - pi^'^+^y^) + log(l - (3^^'^^) 
and therefore to log( ^(xffw )- 

The study of relations between the loop measure and the zeta function (det(/ — 
sP))~^ and more generaly (det(/ — M/P))~^ with / a function on [0, 1] can be done in the 
context of discrete loops. exp(^ based -tWs^'^^VIO) = (det(J — sP))~^ can be viewed 

discretdoops 

as a type of zeta function defined for s G [0 

Primitive (discrete) non trivial based loops are defined as (discrete) based loops which 
cannot be obtained by the concatenation of n > 2 identical based loops. 

The zeta function has an Euler product expansion: if we denote by this discrete 
loop defined by the based discrete loop ^, for ^ = (Ci)---Cfe), A*(C°) — P^lPll----Pe,i ^ 
have: 

(det(/-5P))-i=exp( Yl ^f^'''^^^^^^ n (1 - / ^"^'°V(r))-' 

based primitivediscrete 
discretdoops loops 

4 Geodesic loops 
4.1 Reduction 

Given any finite path uj with starting point Xq, the reduced path uu^ is defined as the 
geodesic arc defined by the endpoint of the lift of uj to 'Xxo- 
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Tree-contour like based loops can be defined as discrete based loops whose lift to the 
universal covering arc still based loops. Each link is followed the same number of times 
in opposite directions (backtracking). The reduced path uo^ can equivalently be obtained 
by removing all tree-contour like based loops imbedded into it. In particular each loop I 
based at xq defines an element in T^o- 




Based loop 
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Reduced based loop 

This procedure is an exemple of loop erasure. In any graph, given a path u;, the loop 
erased path uj^^ is defined by removing successively all based loops imbedded in the path, 
starting from the origin. It produces a self avoiding path (and we see geodesies in 
are self avoiding paths). Hence any non ramified covering defines a specific reduction 
operation by composition of lift, loop erasure, and projection. 




Loop erasure 
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4.2 Geodesic loops and conjugacy classes 

Then, we can consider loops i.e. equivalence classes of based loops under the natural shift. 

Geodesic loops are of particular interest. Note their based loops representatives have to 
be "tailess": If 7 is a geodesic based loop, with I7I = n, the tail of 7 is defined as 7172. ..7? 
if i = sup(j, 7i72...7j = 7„7„_i...7„_j+i). The associated geodesic loop is obtained by 
removing the tail. 

The geodesic loops are clearly in bijection with the set of conjugacy classes of the fun- 
damental group. Indeed, if we fix a base point Xq, a geodesic loops defines the conjugation 
class formed of the elements of r^^, obtained by choosing a base point on the loop and a 
geodesic segment linking it to xq. Any non trivial element of T^o can be obtained in this 
way. 

Given a loop, there is a canonical geodesic loop associated with it. It is obtained by 
removing all tails imbedded in it. It can be done by removing one by one all tail edges 
(i.e. pairs of consecutive inverse oriented edges of the loop). Note that after removal of a 
tail edge, another tail edge cannot disappear, and that new tail edges appear during this 
process. 
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Loop and associated geodesic loop 

A closed geodesic based at Xq is called primitive if it cannot be obtained as the con- 
catenation of several identical closed geodesic, i.e. if it is not a non trivial power in Txo- 
This property is clearly stable under conjugation. Let ^ be corresponding set of primitive 
geodesic loops. They represent conjugacy classes of primitive elements of F (see [37] ). 

4.3 Geodesies and boundary 

Geodesies lines (half-lines) on a graph are defined as paths without backtracking indexed 
by Z (N). 
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Paths and geodesies ean be defined on (X, E) or on a universal cover T and lifted or 

projected on any intermediate covering space. 

Equivalence classes of geodesies half lines for the confluence relation define the bound- 
ary d% of 1. A geodesic half-line on % is can be identified with its origin o and a boundary 
point 6. It projects on a geodesic half-line on X. 

There is a natural cr-field on the boundary generated by cylinder sets Bg defined by 
half geodesies starting with a given oriented edge g. 

Given any point xq in T, one can define a probability measure on the boundary called 
the harmonic measure and dentoted ly^g: ^^^{Bg) is the probability that the lift of the 
P-Markov chain starting at Xq hits g'^ after its last visit to g~. 

Note that F acts on the boundary in such a way that Yi^x) = t^-/{x), for all 7 e T^^. 

Clearly, in the case of a regular graph, as the universal covering is a r— regular tree, 
VxQ^Bg) = ^ (^j._i^d{xo.g-) where d denotes the distance in the tree. When conductances are 
all equal, it can also be computed and is in general distinct from the visibility measure 
from xq, h'^'^''^{Bg), defined as ^ d ^-i ; 2:2... being the points of the geodesic segment 
linking Xq to g~ . 

Note that ^^.g-j S°iy'^^^ [dO] is a shift and F-invariant measure on the set of half-geodesics 
of T. 

It induces a shift and F-invariant probability on half geodesies on X obtained by 
restricting the sum to any fundamental domain and normalizing by \X\. 

4.4 Closed geodesies and associated Zeta function 

Recall that ^ denotes the set of primitive geodesic loops. 
Ihara's zeta function IZ{u) is defined for < < 1 as 

IZ{u) = Y[{1 - vl'^''^)-^ 
7eq3 

It depends only on the graph. 

Note that = E.eq^S?? = E,e<p E^LiPIt)^"^^^) = T^NmU^ where 

denotes the number of tailess geodesic based loops of length m. Indeed, each primitive 
geodesic loop 7 traversed n times still induces ^(7) distinct tailess geodesic based loops. 
Therefore IZ{u) can also be written as exp(E^=2 ^"ilT ) 

Similarly, one can define ITF to be the set of primitive elements of the fundamental 
group F and the F-zeta function to be: 
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TZ{u) = n ~ u^^^Y^ 
7enr 

Note that u '^yz{u) ~ ^mU^ where Lm denotes the number of geodesic based 
loops of length m. TZ{u) can also be written as exp(^^^2 '^"^"' )- Recall that A denotes 
the adjacency matrix of the graph, and D the diagonal matrix whose entries are given by 
the degrees of the vertices. 

Assume now that u < inf(j^). 

Theorem 22 a) '£LmU'^ = {1 - u'^)Tr{[I + {D - I)u'^ - uA]-^) - \V\ 

b) IZ{u) = (1 — u^)~^ det(/ — uA + u'^{D — where x = denotes the Euler number 
\E\ — \X\ of the graph. 

Proof. We adapt the approach of Stark- Terras ([37]) 

It follows from proposition [H] that jVl + ^^L^M'" = {l-u)Tr{G) = {1 -u'^)Tr{[I + 
{D - I)u^ - uA]-^) 

Given a geodesic loop / (possibly empty) and a base point y of I, let Sx,y,i be the sum 
of the coefficients u^'^^\ where 5 varies on all geodesic loops based at x composed with I 
and a tail ending at y. If x = ?/, we have Sy^y^i = u^^^\ Set Sx,i = Sx,y,i- 

Clearly, for any section i of the canonical projection from %x onto X, 

E,,/ Sx,y,i = E,er.-{/} = C-- u)G-'^ - 1 

Therefore, considering first the tailess case, then the case where the tail has length 1, 
and finally decomposing the case where the tail has length at least two according to the 
position of point next to x, we obtain the expression: Sx,y,i = u^^^"* + {dy — 2)^^*^'^+^ + 

= n^(')-«^(')+' + E.(4-i)«'5.,,,/ 

summing in it comes that 

Sx^i = p{iKup(^^ - «^(')+') + E.(4 - iWSx,i 

Then, summing on all geodesic loops / 

EACu - «)^''' - 1) = (1 - «')(E Nmun + - iWiC- - u)G^^^ - 1) 

Therefore, ^ iV„M™ = Tr((/ - u^{D - /)([/ + {D - I)u^ - uA]-'^ - (1 - u'^y^I)) 
= Tr{I - {2u\D -I)- uA)[I + {D - I)u^ - uA]-^ - (1 - u'^)-\I - u\D - I)) 
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= Tr{{2u\D - J) - uA)[I + {D - I)u^ - uA]-^ + (1 - u'^y\u\D - 21)) 
To conclude note that ■^\og{det{I + [D - I)u'^ - uA)) = Tr{{2u{D - 1) - A)[I + [D - 
I)u'^ - uA]-^) and that Tr{u'^{D - 21) = 2u^x- ■ 

Other proofs can be found in the htterature, especially the following one due to Kotani- 
Sunada ([TT]): 

On the line graph, we define a transfer operator Q by = (^^'liz^x}- Then 

IZ{u) = det{{I-uQ)-^ 

Define T and r on A by Ta{x) = J2y {x y}eE '^(^' ra(e) = a(— e). 

Denote S the operator from functions on the graph to A defined by Sf{x,y) = f{y). 

Note that TtS = D,TS = A, and Q = -t + ST. 

Then, for any scalar u, 

(/ - ut){I - uQ) = (1 - n^)/ - (/ - ut)uST 
and 

(/ - uQ){I - ur) = (1 - u^)I - ST{I - ut) (14) 

Therefore T{I - ut){I - uQ) = ((1 - u^)I - uT{I - ut)ST = (/ + u'^iD -I)- uA)T 
and 

T(J - ut){I - uQ){I - ut) = (/ + u\D - J) - uA)T{I - ur) 
Moreover (/ - uQ){I - ut)S = S{{1 - u^)I - uT{I - ut)S) and 

(/ - uQ){I - ut)S = S{{I + u\D -I)- uA)) (15) 

It follows from these two last identities that Im(S') and Ker{T{I — ut) are stable under 
{I -uQ){I -ut). 

Note that S is the dual of —Tt: Indeed, for any function / on vertices and a on 
oriented edges, ^^^^y^^EO a{x,y)S f{x,y) = T^^^^y^^EO a{x,y)f{y) = ^yTTa{y)f{y). 
Therefore, dim(Im(S')) + dim{Ker{TT)) = 2 \E\ 

Note also that dim{Ker{TT)) = dim{Ker{T)) = dim{Ker{T{I — ut)) (as u < 1). 

Moreover, except for a finite set of m's, Im(S') fl Ker(T{I — ut) = {0}. Indeed T(/ — 
ut)S = a — uD which is invertible, except for a finite set of m's. 

Note that ( JT4l) implies that (/ — uQ){I — ut) equals (1 — u"^)! on Ker{T{I — ut) and 
that ([IID implies it equals S{{I + u'^{D - I) - uA))S'''^ on Im(S'). 
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It comes that: 

det((/ - Ut){I - UQ)) = (1 - M2)2|i?|-|X| ^g^^j ^ _ J) _ 

On the other hand, det((/ — ut)) = (1 — m^)'-^', which allows to conclude. 

5 Poisson process of loops 
5.1 Definition 

Still following the idea of [15], which was already implicitly in germ in [39], define, for all 
positive a, the Poissonian ensemble of loops Ca with intensity a/x. We denote by P or 
its distribution. 

Recall it means that for any functional $ on the loop space, vanishing on loops of 
arbitrary small length. 

Of course, Ca includes trivial loops. The periods ti of the trivial loops based at any 
point X form a Poisson process of intensity o.^. It follows (Cf [29] and references therein) 
the following 

Remark 23 The sum of these periods ^ ti and the set of "frequencies" (in decreasing 
order) are independent and follow repectively a T{a) and a Poisson — Dirichlet{0,a) 
distribution. 

Note that by the restriction property, = {/ G CaJ ^ D} is a Poisson process of 
loops with intensity /i^, and that is independent of Ca\C^- 

Remark 24 Note also that these Poisson ensembles can be considered for fixed a or as 
a point process of loops indexed by the "time" a. In that case, Ca is an increasing set 
of loops with stationnary increments. We will denote by CP the associated Poisson point 
process of intensity fi (g) da. It is formed by a countable set of pairs {li,ai) formed by a 
time and a loop 
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We denote by VC^ the set of non trivial discrete loops in Ca- Then, 



with iVi") = N^{C^) = Ei^c^N,{l) and ivg = iV.,,(£a) = Eze£.^x,y(0, when these 
loops are distinct. 

We can associate to a a-finite measure (in fact as we will see, finite when X is 
finite, and more generally if G is trace class) called local time or occupation field 



Then, for any non-negative measure x on X 



and therefore by proposition [17] we have 

Corollary 25 E(e-(^->^)) = det (/ + M^GM^)"" = 



Many calculations follow from this result. 

Note first that E(e-*^" ) = (1 + ^6"^'^)-". Therefore follows a gamma distribution 
r(Q;, G^'"^), with density l{x>o}^r^~ particular, an exponential distribution of 

mean G^'^ for a = 1, as under ^^). When we let a vary as a time parameter, 
we get a family of gamma subordinators, which can be called a "multivariate gamma 
subordinator"[J. 

We check in particular that E(£q, ) = aG^'^ which follows directly from ^{Ix) = G^'^. 

Exercise 26 // = {h}, check that the set of "frequencies" follows a Poisson- 
Dirichlet distribution of parameters (0, a). 

subordinator is an increasing Levy process. See for example reference [1].. 
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(Hint: use the /i-distribution of l^) 
Note also that for a > 1, 



E((l-exp(-^))-i) = C(«). 
More generally, for two points: 

E(e-*^'^e-"^') = ((1 + + 30^'^) - st(G'^'^)2)-" 

X y 

This allows us to compute the joint density of Ca and Ca in terms of Bessel and 
Struve functions. 

We can condition the loop configuration by the set of associated non trivial discrete 
loops by using the restricted cr-field a{VCa) which contains the variables Nx,y. We see 
from [11] and [9] that 



E(e-(^'^)|r'i:«) = JJ( 



K + Xx ' 



The distribution of {N^\x G X} follows easily, from corollary [25] in terms of generating 
functions: 



E(n^^"^^) = det(5.„ + JMi^G.,,. (16) 

X \l \^ Sy 

so that the vector of components A^"'' follows a multivariate negative binomial distribution 
(see for example [S]). 

It follows in particular that Nx"^ follows a negative binomial distribution of parameters 
—a and ^ . Note that for a = 1, Nx^^ + 1 follows a geometric distribution of parameter 



Note finally that in the recurrent case, with the setting and the notations of subsection 
\3.7\ denoting Cae"^ the Poisson process of loops of intensity eajj,^'^'^\ we get that the 
associated occupation field converges towards a constant random variable following a 
Gamma distribution. 
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5.2 Moments and polynomials of the occupation field 

It is easy to check (and well known from the properties of the gamma distributions) that 

^ (a) 

the moments of are related to the factorial moments of Nx ■ 



mc. rm.) ,,,, 

Exercise 27 Denoting C'^ the set of non trivial loops in Ca, define 

k m 
m=l ki+...+k^=k l^^l2...^UeC+ i=l 

Deduce from exerciseMthat E{C~J'^'')vC^) = j^l{N,>k}{N^"^ -k + l)...(A^i"^ - l)Nt^ 



It is well known that Laguerre polynomials L^^ with generating function 



j:t^Lt'\u) 



e 



(i-ty 



are orthogonal for the T{a) distribution. They have mean zero and variance ^-^^|p^- Hence 
if we set = G^'^and P^''^{x) = {—a)^ L^i^~^\^) , the random variables P^''^'"(£a ) are 
orthogonal with mean and variance cr^^ ^^"k\^^ ? fo^' ^ > 0. 

Note that Pi'"''{Ca ) = Ca — QcTx = Ca — ]E(£q, ). It will he denoted Ca ■ 

ut 

Moreover, we have Eot'Pn^) = E(-^t)'4""'^(^) = ^ 
Note that by corollary! 



E' 



1 + or^t)" (1 + ays) 



(1 + (7a;t)°(l + CTj^s)" l + (7xt l + ayS 1 + axtl + CTyS 

= (1 - st(G^'«')2)-". 

Therefore, we get, by developping in entire series in (s,t) and identifying the coefficients: 
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Let us stress the fact that G^'^ and G^'^ do not appear on the right hand side of 
this formula. This is quite important from the renormahsation point of view, as we will 
consider in the last section the two dimensional Brownian motion for which the Green 
function diverges on the diagonal. 

More generally one can prove similar formulas for products of higher order. 

Note that since G^M^ is a contraction, from determinant expansions given in [IQ] and 
|41] . we have 

oo 

det(/ + M^GM^)-" = 1 + J2 Xn■■■X^,Pera{G.,,^^, l<l,m<k) (18) 

k=l 

and then, from corollary [251 it follows that: 

H{^c.,x) ) = X^Xn---Xifc^er„(Gi„i„, l<l,m<k) 

Here the a-permanent Per^ is defined as ^o-g^^ '^'"*''^^^n,ia(i)---^jfc,V(fc) with m{a) denoting 
the number of cycles in a. 

Note that an explicit form for the multivariate negative binomial distribution, and 
therefore, a series expansion for the density of the multivariate gamma distribution,follows 
directly (see [H]) from this determinant expansion. 

It is actually not difficult to give a direct proof of this result. Thus, the Poisson process 
of loops provides a natural probabilistic proof and interpretation of this combinatorial 
identity (see [S] for an historical view of the subject). 

We can show in fact that: 

Proposition 28 For any (xi, ...Xk) in , E(£q: ^ ---Ca = PeVaiG^^'^"" , 1 <l,m < k) 

Proof. The cycles of the permutations in the expression of Pera are associated 
with point configurations on loops. We obtain the result by summing the contributions 
of all possible partitions of the points ii...ik into a finite set of distinct loops. We can 
then decompose again the expression according to ordering of points on each loop. We 
can conclude by using the formula /^(/^^'•••'^'") = G^^'^'^G^'^'^^ ...G^"^'^^ and the following 
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property of Poisson measures (Cf formula 3-13 in [12]): For any system of non negative 
or yU-integrable loop functionals Fj 



Remark 29 We can actually, in the special case ii = i2 = ■■■ = ik = x, check this 
formula in in a different way. From the moments of the Gamma distribution, we have that 
E((£q )") = {G^'^)^a{a+l)...{a+n—l) and the a-permanent can he written 'Y^^d{n,k)a^ 
where the coefficients d{n,k) are the numbers of n— permutations with k cycles (Stirling 
numbers of the first kind). One checks that d{n + 1, A;) = nd{n, k) + d{n, k — 1). 

Proposition 30 Given any bounded functionals $ on loops configurations and F on 
loops, we have: 

EiJ2 nm^a)) = I U {l}))aF{l)ii{dl) 

Proof. This is proved by considering first for the functionals of the form 

Tl,l-^Jkl2 f^iqeCc 111 ^ji^j)) (with Gj bounded and /i-integrable) which span an algebra sep- 
arating distinct configurations and applying formula (fT9|) : Then, the common value of 
both members is a" KFGj) Ui^^ fi{Gi) + a^+^fi{F) HI l^{Gj) ■ 

The above proposition applied to F{1) = l^,Nx"y and propositions [T2] and [T8] yield the 
following: 

Corollary 31 EmC^)TJ) = aj E^C^Ui^jW ij{d^) =aj EmC^U{^}))iJ^'^{d^) 
and 

E($(£„)ivS) = a jEi^C^ U {7}))iV,,,(7)/i(d7) = aG^^yjEi^C^ U {tDK-^^^t) 
ifxj^y. 

Let be the set of permutations of k elements without fixed point. They correspond 
to configurations without isolated points. 

Set Per°(G*''*'", 1 < /,m < A;) = E^e^o a"'^''^G'^''-w ...G'^''^") . Then an easy calcula- 
tion shows that: 
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Corollary 32 ¥.{C^\..C^'') = PerO(G^''^-, l<l,m<k) 
Proof. Indeed, the expectation can be written 

P<k /c{i,...fe},|/|=p ie/= 

and 

Per«(G'^-^^ a,beI) = J2 H G^'^ Per1{G'^''\ a.beJ). 

^C/ j(^I\J 

' — i\ ' — i]^ „ 

Then, expressing E(£q, ...£q, ) in terms of Perm's, we see that if J C {1, ...k), \ J\ < k, the 

coefficient of Per° a, 6 G J) is ^ JJ G'^''^ which vanishes as (-l)"'^' = 

I,IDJ jeJ" 
= (_i)l^l(_i)IAJ| and E/Dj(-1)''^^' = (1 - 1)'"'^' = 0. ■ 

Set Q'^''^{u) = P^''^{u + aa) so that P^'"{Ca ) = Q'^''^{Ca )• This quantity will be 
called the n-th renormalized self intersection local time or the n-th renormalized power 
of the occupation field and denoted 

From the recurrence relation of Laguerre polynomials 

nL^^-'\u) = {-u + 2n + a~ 2)4^/^ - {n + a - 2)L^^:^\ 

we get that 

(^) = - Mn - l))QZi{n) - cy\a + n - 2)Qr^^[u) 

In particular Q'^'" [u] = — 2au — aa^). 
We have also, from (JTTl) 

E(Qr^(Z;^),Qrn2;')) = 6,,iG^'yr^^^±^^^^^^±^^ (20) 

The comparison of the identity fl20l) and corollary [32] yields a combinatorial result 
which will be fundamental in the renormalizing procedure presented in the last section. 

The identity (120|) can be considered as a polynomial identity in the variables cTj;, ay 
and G'^'y. 

Set <5fc''^^(M) = ^m=o Qm^u^f^t'"^ deuote Nn,m,r,p the number of ordered configu- 
rations of n black points and m red points on r non trivial oriented cycles, such that only 
2p links are between red and black points. We have first by corollary [321 
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r p<'mf{m,n) 

and therefore 

J2 Yl «'C'9:''^n,m,.,p = unless p = I = k. (21) 

r p<m,<k p<n<l 

r a,k a,k^r a(a + 1) ■ . . (q + fc - 1 ) 

^fe ^fe ^fc,fc,r,fe = (22) 

r 

Note that one can check directly that g^'^ = and A^fc,fc,i,fc = ^'(^ " l)'; Nk,k,k,k = ^! 
which confirms the identity (122!) above. 

5.3 Hitting probabilities 

Denote by [H^]y = F^{xtp = y) be the hitting distribution of F by the Markov chain 
starting at x {H^ is called the balayage or Poisson kernel). Set D = F'^ and denote by , 
= P\^^j^^ = [{I - P^)]-i and = [(Ma - C)\dxd]~^ the energy, the transition 
matrix, the potential and the Green function of the Markov chain killed at the hitting 
time of F. 

Denote by the law of the killed Markov chain starting at x. 
Recall that 

[H^Ty = + E.eD[iP''n.Py = M^=y} + E.eDl^'r.Py ■ Moreover we have 

by the strong Markov property, V = + H^V and therefore G = G^ + H^G. (Here 
we extend and G^ to X x X by adding zero entries outside D x D). 

As G and G^ are symmetric, we have [H^G]y = [H^G]^ so that for any measure u, 
H^{Gu) =G{uH^). 

Therefore we see that for any function / and measure z/, e{H^f, G^v) = e{H^f, Gv) — 
e{H^f, H^Gu) = {H^f, v) - e{H^f, G{H^u)) = as {H^f = . 
Equivalently, we have the following: 

Proposition 33 For any g vanishing on F, e{H^f,g) = so that I — is the e- 
orthogonal projection on the space of functions supported in D. 

H^l is called the capacitary potential of F and e{H^l, H^l) the capacity of F. 
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Note that these results extend without difficulty to the recurrent case. In particular, 
for any measure v supported in G^v = G{v — vH^) and e{H^ f, G^u) = for all /. 
For further developments see for exemple ( [18]) and its references. 

The restriction property holds for as it holds for /x. The set of loops inside D 
is associated with /i^ and is independent of Ca — i^a- Therefore, we see from corollary 
[2S]that 

-(r.-H5,x). ^ det(GJdet(G^) 
^ ^ ^ det(G') det(G^)^ " 

From the support of the the Gamma distribution, we see that > 0) = oo. But 

this is clearly due to trivial loops as it can be seen directly from the definition of /i that 
in this simple framework they cover the whole space X. 

Note however that 

lxff{F) > 0,p > 1) = fi{p > 1) - = 0,p > 1) = ^(p > 1) - > 1) 

= -log( ^^^(^-^^ ) = -log( ^^^^^"^ ). 

It follows that the probability that no non trivial loop (i.e. a loop which is not reduced 
to a point) in intersects F equals 

eM-MMD > i,m > 0})) = 

Recall that for any {n + p,n + p) invertible matrix A, denoting Cj the canonical basis, 

det(y4"^) det{Aij, I <i,j <n) = det(A"^) det(Aei, ...Ae„, e„+i, ...e„+p) 

= det(ei, ...e„, A ^e„_|_i,...A ^e„+p) 
= det{{A-^)k,i,n< kj<n + p). 

In particular, det(G^) = ^'^^^^ , so we have the 

det(G|_FxF) 

Proposition 34 The probability that no non-trivial loop in Ca intersects F equals 

[TT A, det (G)]-". 

'■J-J- FxF ' 

Moreover E{e<^-^^^)) = ( d^t^xH^xK . 
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In particular, it follows that the probability that no non-trivial loop in visits x 
equals qx,x )'^ which is also a consequence of the fact that follows a negative binomial 
distribution of parameters —a and j^^^ 

Also, if Fi and F2 are disjoint, 

li{T{Fi)T{F2) > 0) = > 0,p > 1) + ^x{T{F2) > 0,p > 1) - ii{T{Fi U F2) > 0,p > 1) 

det(G') det(G'^i^^^) 
~ det(G'^i)det(G'^2) >■ 

Therefore the probability that no non-trivial loop in Ca intersects Fi and F2 equals 

.X ttO s ,det(G)det(G^in^2) 

It follows that the probability no non-trivial loop in Ca visits two distinct points x and y 
equals ( '^ ' qI'I^^^ '^^ )" and in particular 1 — ^,xQy,y if a = 1. 

Exercise 35 Generalize this formula to n disjoint sets: 

det(G)n,<<det(Ga™>) 



Note this yields an interesting determinant product inequality 



6 The Gaussian free field 

6.1 Dynkin's Isomorphism 

By a well known calculation, if X is finite, for any % G M^, 



(27r) 1^1/2 y y ^g^^^^ 

and 



ydet(MA - C) /•^^^^^_i<,2,^>^_i,(,)^ ^„u_.^Na:,2/ /^^K^xy 
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This can be easily reformulated by introducing on an independent probability space 
the Gaussian field (f) defined by the covariance E<^(0'^0^) = G^'^ (this reformulation cannot 
be dispensed when X becomes infinite) 

So we have E<^(e-5«^''X>) = det(/ + GMJ-^ = ^det(G'^G-i) and 
E0((0'^0^e~2'^''^^'^^) = {G^^^y ^J(\.ei{G^G^^) Then since sums of exponentials of the 
form e~^^ ''^^ are dense in continuous functions on the following holds: 

Theorem 36 a) The fields Ci and have the same distribution. 

b) E0((0^0^F(i0^)) = jE(F(£i +^))fi^'y{d'j) for any bounded functional F of a non 
negative field. 

Remarks: 

a) This is a version of Dynkin's isomorphism (Cf [7]). It can be extended to non- 
symmetric generators (Cf |20j). 

b) By corollary [3ll if Cx^y 7^ 0, b) implies that 



E,((0>^F(^0^)) = 7^ / E(F(£j))ivi|) 



c) An analogous result can be given when a is any positive half integer, by using real 
vector valued Gaussian field, or equivalently complex fields for integral values of a (in 
particular a = 1): If = (0i, 02---0fc) are k independent copies of the free field, the 

fields Ck and | = ^Yli'P'j have the same law and (p^^ F{^\\(p\f)) = 

k /E(F(£l +7))/i^'2^(rf7) 

d) Note it implies immediately that the process 0^ is infinitely divisible. See |8] and 
its references for a converse and earlier proofs of this last fact. 

Remark HT] and theorem [36] suggest the following: 

Exercise 37 Show any bounded functional F of a non negative field, if Xi are 2k points: 

E(F(0^)n<^")=/E(^£+ET;-)) E n/^'^'^^(^^^-) 

1 pairings pairs 

where Ylpairings '^^^'^^ that the k pairs yj,Zj are formed with all the 2k points Xi, in 
all ^IIp possible ways. 
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Hint: As in the proof of theorem we take F of the form e 2< 'X> and prove the 
identities from formula [29] by perturbation and differentiation. 

Exercise 38 For any f in the Dirichlet space M of functions of finite energy (i.e. all 
functions if X is finite ), the law of f + (p is absolutely continuous with respect to the law 
of (j), with density exp(< —Lf,(j) — ^e(/)) 

Exercise 39 a) Using proposition show (it was observed by Nelson in the context 
of the classical (or Brownian) free field) that the Gaussian field cj) is Markovian: Given 
any subset F of X, denote Hp the Gaussian space spanned by {(p^.y G F}. Then, for 
X E D = F^, the projection of (j)^ on Tip (i-e. the conditional expectation of (j)^ given 
a{cpy,yeF))is Y.y^p[H^YyCt>y . 

b)Moreover, show that (p^ = (p — H^cp is the Gaussian field associated with the process 
killed at the exit of D. 

6.2 Fock spaces and Wick product 

The Gaussian space Ti spanned by {(p^,x G X} is isomorphic to the Dirichlet space EI 
by the linear map mapping 0^ on G^'-. This isomorphism extends into an isomorphism 
between the space of square integrable functional of the Gaussian fields and the symmet- 
ric Fock space r{M) = ©H*^" obtained as the closure of the sum of all symmetric tensor 
powers of H. This is Bose second quantization: See for a complete description 

of this isomorphism. 

We have seen in theorem [36] that functionals of Ci can be represented in this 
space of Gaussian functionals. In order to prepare the extension of this representa- 
tion to the more difficult framework of continuous spaces (which can often be viewed 
as scaling limits of discrete spaces), including especially the planar Brownian motion 
considered in [15], we shall introduce the renormalized (or Wick) powers of (p. We set 
: ((/)^)" := {G^'^)2Hn{(p^/vG^) where Hn in the n-th Hermite polynomial (characterized 
by Yl \iiHn{u) = e*"~~). It is the inverse image of the n-th tensor power of G^'' in the 
Fock space. 

Setting as before = G^'^, from the relation between Hermite polynomials and 

_ 1 

Laguerre polynomials L„ ^ , 
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H2n{x) = (-2)"n!L;^(y) 

it follows that: 

More generally, ii (j)i,(j)2---4>k are k independent copies of the free field, we can define 
• 11^=1 • = n^=i • •• Then it follows that: 

1 ni + ..+nfc=n j=l 

From the generating function of the polynomials Pn ' , 

k 

1 

Therefore, 

rf^"(^) = ^ (23) 

Note that : 0^ := 0^ — 0" These variables are orthogonal in L^. Let = P — a be 
the centered occupation field. Note that an equivalent formulation of theorem [36] is that 
the fields | : Yl\ 0j • and Ck have the same law. 

Let us now consider the relation of higher Wick powers with self intersection local 
times. 

Recall that the renormalized n-th self intersections field C^'"^ = P^'"{Ca ) = Q^^i^^a ) 
have been defined by orthonormalization in of the powers of the occupation time. 
Then comes the 

Proposition 40 The fields C'C and : '/'j)" ■ have the same law. 

This follows directly from fl23l) . 



^ ni\...nk\ 1 ^ J' 

ni+..+rn.=n j=l 
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Remark 41 As a consequence, it can be shown that: 

r 

where Ski^k2,...kj is the set of permutations a of k = ^ kj such that 

c^({Er^ + !>••• + ^j}) ^ {Er^ + 1> ••• Zir^ ^« + ^j} empty for all j. 

The identity follows from Wick's theorem when a is a half integer, then extends to all 
a since both members are polynomials in a. The condition on a indicates that no pairing 
is allowed inside the same Wick power. 

7 Energy variation and representations 

The loop measure /i depends on the energy e which is defined by the free parameters 
C,K. It will sometimes be denoted /ig. We shall denote Ze the determinant det(G') = 
det(M, - C)-\ Then > 1) = log(Ze) + Elog(A.). 
is called the partition function of Ca- 
The following result is suggested by an analogy with quantum field theory (Cf [115]). 

Proposition 42 = 

n) IfC.,y >0,J^ = -T.,yf^ wtth T^,yil) = i?^+Ty) - ^(l) - ^(l). 

Proof. Recall that by formula ©: fi*{p = l,C = x,t e dt) = e~^^*y and = 

k, = x„ f, G dt,) = I u.,y c^r n. a-^^ n^ezM h.^-'^^'^dt. 

Moreover we have C^^y = Cy^^ = X^Py and = + Y.y Cx,y 
The two formulas follow by elementary calculation. ■ 
Recall that = G^'^ and f^{N^,y) = C^'yC^^y. 
So we have = G^'"^ + G^'-f - 2G^'^. 

Then, the above proposition allows us to compute all moments of T and / relative to He 
(they could be called Schwinger functions). 

Exercise 43 Use the proposition above to show that: 

j Tjyfi{di) = {G^'yf 
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y"r,T,,,(/)/x(c?/) = (G^'^-G^'^)2 

and 

j Tx^y{l)Tu^vi^) fJ'{dl) = {Gx,U + Gy^v — Gx,V — Gy^uY 

The calculations are done noticing that for any invertible matrix function M(s), 
f^M{s)-^ = -M{s)-^M'{s)M{s)-\ The formula is applied to M = Ma - C and s = 
or Gx,y 

We can apply these results to the Poissonnian loop ensembles, to get the following 
Corollary 44 For any bounded functional $ on loop configurations 

ii) IfGx,y > 0, ^E($(/:„)) = E(tS<I>(/:,)) with r^'j = Eiec^T^^yi^) = i^a + 

^D-tB'S '"'"^ ^^"^y = ^^"^^ ~ ^^^^^i) = - "(G"'" + Gy^y - 2G-^y). 

The proof is easily performed, taking first $ of the form X^ii^Za i-iq<^Cc Hi ^i(^i))- 
The proposition |12] is in fact the infinitesimal form of the following formula. 

Proposition 45 Consider another energy form e' defined on the same graph. Then we 
have the following identity: 

Consequently 

fj^eiie^ ^''^y - 1)) = log( — ) (24) 

Proof. The first formula is a straightforward consequence of ([7j). The proof of ^IM 
goes by evaluating separately the contribution of trivial loops, which equals ^a.log(^). 
Indeed, 

^^((eE^^.i°^(if)-M-A.)P_^) = ^,,{p>l)-^,^{p>l)+ ;.,(l^^^,j(eE(Ai-A.)p_i)). 
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The difference of the first two terms equals \og{Ze')+Y^ log(A'^) — (log(Ze)— ^ log(A^)). 
The last term equals J^x /o°°(^ ~ l)^dt which can be computed as before: 

/ie(l{p=i}(e^(^-^^)'^ - 1)) = -J2^og& (25) 



Remark 46 ( h-transforms) Note that if C'^ y = h^h^Cx^y and k'^ = —h^{LhyXx for 
some positive function h on E such that Lh < 0, as X' = h'^X and [P']y = -^Pyhy , we 
have [Gr'y = §j^and%^ = ^. 

Remark 47 Note also that [^]^ = E(e~5[<='"'=K9^)^^ 2s the Gaussian free field associ- 
ated with e. 



Integrating out the holding times, formula fl2^ can be written equivalently: 

f^eiUl^f-'Uiv^'''''' - 1) = (26) 



and therefore 



/ N Gxy X Cxy Zf, 

{x,y) X {x,y) 



Note also that YlixJij^r-^ = U^x.JlT^]''-^^''-^ ■ 



Remark 48 These ^ determine, when e' varies with ^ < 1 and y = 1, the Laplace 
transform of the distribution of the traversal numbers of non oriented links Nx,y + Ny^x ■ 

Other variables of interest on the loop space are associated with elements of the space 
of odd real valued functions uj on oriented links : uj^'"^ = —uj"^'^. Let us mention a few 
elementary results. 
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The operator [-P*-'^-']^ = Py exp{iuj^'y) is also self adjoint in L'^{\)- The associated loop 
variable can be written yUJ^'^Nx^y{l). We will denote it f^u. This notation will be 
used even when u is not odd. Note it is invariant if u is replaced hj u + dg for some 

g. Set [G'^'^^]^'^ = -. By an argument similar to the one given above for the 

occupation field, we have: 

n,.(e^^'- - 1) = exp(t(P(-) - - exp(t(P - J))^. 

Integrating in t after expanding, we get from the definition of /i: 

j {e^f'- - l)df,il) = ^[Tri{P^-Y) - Tr{{Pn. 

Hence /(e^/^'^ - l)dfi{l) = log[det(-L(/ - pH)-i] = log(det(G(^)G-i)) 

We can now extend the previous formulas (12^ and (ITTj) to obtain, setting det(G'^'^)) = 



■'e,iL> 



gEA^.,.iog(^)-E(AL-A.)r.+./,. _ i-^^^^^i-^ ^ log(%^) (28) 

Ze 

and ^ 

^^-Q[^gic...,|7vWg-E{Al-A.)2;^) = (29) 
^x,v Z^ 

x,y 

Remark 49 The a-th power of a complex number is a priori not univoquely defined as a 
complex number. But log[det(/ — P^^^)] and therefore log(-Ze,(^) o'^e well defined as P^^^ 
is a contraction. Then Z^^ is taken to be exp(a log(2e,a;))- 

To simplify the notations slightly, one could consider more general energy forms with 
complex valued conductances so that the discrete one form is included in e'. But it is 
more interesting to generalize the notion of perturbation of P into P^"^^ as follows: 

Definition 50 A unitary representation of the graph (X, E) is a family of unitary ma- 
trices [[/^'^], with common rank djj, indexed by E'-^ , such that [U^'^] = [[/^'^]~^. 
We set P(^) = P®U (more explicitly [P^^^xi = 

Similarly, we can define C^^^ = ^P^^\ V^^^ = (/ - P^^Y\ G^^^ = 

One forms define one- dimensional representations. The sum and tensor product of 

two unitary representations U and V are unitary representations are defined as usual, 

and their ranks are respectively du + dy and dijdy- 
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Definition 51 Given any based loop I, set tjj{1) = -^Tr(WU^^'^^+^) if p{l) > 2 and the 
associated discrete based loop is {C,i,C,2--,(,p), and tu{1) = 1 if p{l) = 1. 
For any set of loops C, we set tu{C) = Yliec '^ui})- 



Remark 52 a)\Tu{l)\ < 1 

b) Tu is obviously a functional of the discrete loop ^ 

c) tu{1) = 1 if C,° is tree-like. In particular it is always the case when the graph is a 
tree. 

d) If U and V are two unitary representations of the graph, tjj^v = ^i/ + Ty and 
Turn = I'U'TV- 

From b) and c) above it is easy to get the first part of 

Proposition 53 i)The trace tu{1) depends only on the canonical geodesic loop associated 
with the loop ^° , i. e. of the conjugacy class of the fundamental group associated with the 
element o/F^i defined by the based loop ^. 

ii)The variables tu{1) determine, as U , varies, the geodesic loop associated with I. 

Proof. The second assertion follows from the fact that traces of unitary represen- 
tations separate the conjugacy classes of finite groups (Cf ) and from the so-called CS- 
property satisfied by free groups (Cf t38]|): Given two elements belonging to different 
conjugacy classes, there exists a finite quotient of the group in which they are not conju- 
gate. ■ 

Again, by an argument similar to the one given for the occupation field, we have: 
K^ru - 1) = i E£i exp(t(p(^) - 7))^;: - exp(t(P - 
Integrating in t after expanding, we get from the definition of /i: 

oo 

k=i ^ 

We can extend P into a matrix P^^^'u^ = p ^ J^^ indexed by X x {1, du} by taking its 
tensor product with the identity on M.'^^ . 
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Then: 

^ k=l 

Hence, as in the case of the occupation field jijjjil) — l)djji{l) = ^ log(det(l^*^'^)) 
= log(det(G(^))) - log(det(G')) as det(G (S) hu) = du det(G) 

Then, denoting 2^e,c/ the ^-th power of the determinant of the (|X| du, \X\ du) matrix 
G(^) (well defined by remark l49l) . the formulas (|28l) and (!29l) extend easily to give the 
following 



C / ^ 

Proposition 54 a;/(e^^^'^'°'^^^^"^^^-~^^^^V[/(/) - l)fie{dl) = log'^^' 



e'.U ~ 



Let us now introduce a new 



Definition 55 We say that sets Aj of non-trivial loops are equivalent when the associated 
occupation fields are equal and when the total traversal numbers YIi^a ^x,y{l) <2?^e equal 
for all oriented edges {x,y). Equivalence classes will be called loop networks on the graph. 
We denote A the loop network defined by A. 

Similarly, a set L of non-trivial discrete loops defines a discrete network characterized 
by the total traversal numbers. 



Note that these expectations determine the distribution of the network Ca defined by 

C f ^ 

the loop ensemble £„. We will denote 5^''^'''^ (/) the variables e^^'''''*''^'"^*"'^^'!''' '^^'^'^ \x)ix+iSit^ 
and B^'^''^{Ca) the variables 



n B-^'-{i) = 



More generally, we can define B^''''^{1) and B''''''^{Ca) in a similar way as B"'^ '^(/) and 
B^''^'''^{Ca), using Tu{l) instead of e*-^;"^. Note that for each fixed e, linear combinations 
of the variables B^''^''^(Ca) form an algebra as B^'^'^''^^B^''^2'^^ = B^'^'^^^''^'^'^®^'^ , and in 

particular, ^e^e'i.'^i^e.e^.c^a = ^e,e;+e^,c^i+a;2_ 
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Remark 56 Note that the expectations of the variables B^''^' (Ca) determine the dis- 
tribution of the network La defined by the loop ensemble Ca- The expectations of the 
variables B^'^' [Co) determine in addition the law of the set of geodesic loops associated 
with Ha- 

To work with /i, we will rather consider linear combinations of the form ^ \i(B^'^i'^^ — 
1), with ^ Aj = 0, which form also an algebra. 



Remark 57 Formulas and [29\) apply to the calculation of loop indices: If we have 
for example a simple random walk on an oriented planar graph, and if z' is a point of the 
dual graph X' , u^^'^ can be chosen such that for any loop I, ji^^^''^ is the winding number 

of the loop around a given point z' of the dual graph X' . Then e^^ ^'^'^^ ■''^ is a spin 
system of interest. We then get for example that 



Jl ^TT Jq 

and hence 

P( J2 I fuj^^'^l = 0) = e5^/o^'°s('i^*(«''"""'''^^^"'))'^" 

Conditional distributions of the occupation field with respect to values of the winding num- 
ber can also be obtained. 



8 Decompositions 

Note first that with the energy e, we can associate a rescaled Markov chain Xt in which 
holding times at any point x are exponential times of parameters A^;: Xt = Xrt with 
Tt = inf(s, -^du = t). For the rescaled Markov chain, local times coincide with 
the time spent in a point and the duality measure is simply the counting measure. The 
Markov loops can be rescaled as well and we did it in fact already when we introduced 
pointed loops. More generally we may introduce different holding time parameters but it 
would essentially be useless as the random variables we are interested in are intrinsic, i.e. 
depend only on e. 
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li D G X and we set F = D"^, the orthogonal decomposition of the energy e(/, /) = 
e(/) into e^(/ — f) + e{H^ f) (see proposition [331) leads to the decomposition of the 
Gaussian free field mentioned above and also to a decomposition of the rescaled Markov 
chain into the rescaled Markov chain killed at the exit of D and the trace of the rescaled 
Markov chain on F, i.e. xl^^ = Xgp, with S'f = inf(s, lF(xu)du = t). 

Proposition 58 The trace of the rescaled Markov chain on F is the rescaled Markov 
chain defined by the energy functional e^^^{f) = e{H^f) , for which 

^x^^ — ^x,y + Cx,aCb,y[G^]"'''' 
a,beD 

^ = - E c.MgT' 

a,beD 

and 

Ze = Zf,DZ^{F} 

Proof. For the second assertion, note first that for any y ^ F, 
Moreover, e{H^f) = e(/, H^f) and therefore 

aeD 

where = T.afieD Pa[G^T'''Cb,x = 'EaeD Pa[H^]x is the probability that the Markov 
chain starting at x will return to x after an excursion in D. 
Then for distinct x and y in F, 

Ci? = -e^^Hkx},l{y}) = -e(l|.|,i^^l|,}) 

= Cx,y + CxAH^]y = Cx,y + Cx,aCb,y[G^]""'^ . 

a a,b£D 

Note that the graph defined on F by the non-vanishing conductances Gx,y has in 
general more edges than the restiction to F of the original graph. 
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For the third assertion, note also that G^^^ is the restriction of G to F as for all x,y G 
F, e^^^{G5y\FA{x}) = e{G5y, [H^l^^}]) = ^{x=y}- Hence the determinant decomposition 
already used in section [5^3] yields the final formula. The cases where F has one point was 
already treated in section I5.3[ 

For the first assertion note the transition matrix [-P''^^-'']^ can be computed directly and 
equals 

Py+ Ea.feeD ^^""^"^Ife^i = Py+ Ea.feeD [G'^''^"^]"'^^,^ • It cau be decomposed 
according to whether the jump to y occurs from x or from D and the number of excursions 
from X to x: 

oo 

k=0 a,beD a,beD 
oo 

= T.(T.Pa[G''r''c,,,nPy + J2 p:[G''r''c,,y). 

k=0 a,beD a,beD 

The expansion of in geometric series yields exactly the same result. 

Finally, remark that the holding times of x] at any point x E F are sums of a 
random number of independent holding times of Xt- This random integer counts the 
excursions from x to a; performed by the chain xt during the holding time of xl . It 
follows a geometric distribution of parameter 1 —pi^\ Therefore, —rpr = ^—tft- is the 

expectation of the holding times of xl at x. ■ 

If X is carried by D and if we set = e + || denote [e^^j'^^'^J' by e^^'^^ we 

have 

cis^^ = c.,y + E Cx,aa,y[G^r^ = J2 p^Kr^'c,,^ 

a,b a^baD 

and A™ = A.(l-p™). 

More generally, if is such that = C on F x F, and A = A* on F we have: 

G*P = Cx,y + E ClCl[G*^r, V*J'' = E Pt[G*^rcl 

a,b a,b(iD 

and Af^^ = A,(l-pf^^). 

A loop in X which hits F can be decomposed into a loop l^^^ in F and its excursions 
in D which may come back to their starting point. Let /x^* denote the bridge measure 
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(with mass [G^]"''^) associated with e^. 
Set 

1 oo ^ 

^x,y a,beD n=l '^a;]?a; a,6eZ3 

and = + T^nM^^P^n 

Note that;>f (1) = = (1) = 1. 

A loop I hitting F can be decomposed into its restriction l^^^ — {^i,Ti) in F (possibly 
a one point loop), a family of excursions 7^i,5j_|_i attached to the jumps of l^^^ and systems 
of i.i.d. excursions (7^., h < n^.) attached to the points of Note the set of excursions 

can be empty. 

Let be the restriction of /i to loops in contained in D. It is the loop measure 
associated to the process killed at the exit of D. We get a decomposition of /x — /^^ in 
terms of the loop measure /x^^^ defined on loops of F by the trace of the Markov chain 
on F, probability measures u^y on excursions in D indexed by pairs of points in F and 
measures on excursions in D indexed by points of F. Moreover, the integers follow 
a Poisson distribution of parameter ^\f^Ti (the total holding time in before another 
point of F is visited) and the conditional distribution of the rescaled holding times in 
before each excursion 7|. is the distribution /5„^.,t* of the increments of a uniform sample 
of n^. points in [0 Tj] put in increasing order. We denote these holding times by Ti^h and 
set/ = A(/i^},(76:,e.^J,(n^„7tr,.)). 

Then /x — fi^ is the image measure by A of 

k 

Note that for x, y belonging to F, the bridge measure ii^'^ can be decomposed in the 
same way, with the same excursion measures. 

The one point case 

If F is reduced to a point x, the decomposition is simpler, l^^^ is a one point loop 
with r = — and the number of excursions (all independent with the same distribution 

pf^ ) follows a Poisson distribution of parameter }^^t — X^l^ . 
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The law can of course be decomposed in the same way, with the same conditional 
distribution given l^. Recall it follows an exponential distribution with mean G^'^. 

follows a r(a,G^'^) distribution, in particular an exponential distribution with 
mean G^'^ for a = 1. Moreover, the union of the excursions of all loops of outside 
X has obviously the same Poissonian conditional distribution, given = s than fi and 
The set of excursions outside x defined by the ^^-distributed bridge and by Ci 
are therefore identically distributed, as the total holding time in x. 

Remark 59 Note finally that by exercise [2R the distribution ofCi/C\^^ can be recovered 
from a unique sample of by splitting the bridge according to an independent sample 
Ui of Poisson — Dirichlet{0,a) , more precisely, by splitting the bridge (in fact a based 
loop) I into based subloops l\[ai,ai+i], with cTj = inf(s, j- l{x}{h)ds = Uil^). 

Conversely, a sample of the bridge could be recovered from a sample of the loop set 
Ci/C^^ by concatenation in random order. This random ordering can be defined by 
taking a projective limit of the randomly ordered finite subset of loops {/j,n} defined by 
assuming for example that l^^ > ^. 

Coming back to the general case, the Poisson process cP = G C^} has 

intensity fi^^^ and is independent of C^. 

r p\ 

Note that Ca is the restriction of Ca to F. 
If X is a measure carried by we have: 



E(e-(S;'«>|£m) =E(e-(^*>)( H l/e~'™'''f„(<'7)]"'''"'""' 

n 



x,yeF 

r^ia. TT r^^^^r 'ri^hTT r.{^.x}_,W 



x,y(^F ^x,y 



r p\ 

(recall that Ca is the restriction of to F). Also, if we condition on the set of discrete 
loops VCl[^ 

x,y&F ^x^y x&F '^^ 
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where the last exponent A^^, + 1 is obtained by taking into account the loops which have 
a trivial trace on F (see formula fl25|) ). 

More generally we can show in the same way the following 

Proposition 60 // C* = C on F x F, and \ = on F , we denote B'^'^* the multi- 
plicative functional ||[-^;^]^-^!'e~^-6-f' "'^"^ 

Then, 



C'x,y 



and 



^#{F} 

x,y£F ^^^y x£F 
^ ^#{F} 

x,y£F ^^>y xeF 

These decomposition and conditional expectation formulas extend to include a current 
u in C*. Note that if u is closed (i.e. vanish on every loop) in D, one can define such 
that [Ce*'"]{^> = C^^^e^'^^. Then 

The previous proposition implies the following Markov property: 

Remark 61 If D = DiU D2 with Di and D2 strongly disconnected, (i.e. such that for 
any {x, y, z) E Di x D2 x F, C^^y and C^^zCy^z vanish), the restrictions of the network Ca 
to DiU F and D2U F are independent conditionally on the restriction of Ca to F. 

Proof. It follows from the fact that as Di and D2 are disconnected, any excursion 
measure u^y or from F into D = DiU D2 is an excursion measure either in Di or in 
D2. m 



Branching processes with immigration An interesting example can be given after 
extending slightly the scope of the theory to countable transient symmetric Markov chains: 
We can take X = N — {0}, Cn,n+i = 1 for all n > 1 and ki = 1. P is the transfer matrix 
of the simple symmetric random walk killed at 0. 
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Then we can apply the previous considerations to check that is a branching process 
with immigration. 

The immigration at level n comes from the loops whose infimum is n and the branching 
from the excursions of the loops existing at level n to level n + 1. Set F„ = {1, 2...n} and 

From the calculations of conditional expectations made above, we get that for any 
positive parameter 7, 

Prom this formula, it is clear that is a branching Markov chain. To be more precise, 
note that for any n,m > 0, — 2(n A m) and A„ = 2 and that G^'^ — 1. Moreover, 
Gl^l = G^'^'-G^'^-fG^'g^ so that G^f^ = i^^and for any n > 0, the restriction of the Markov 
chain to D„ is isomorphic to the original Markov chain. Then it comes that for all n, 
p{^"> = 1, Ai^"^ = 1, and Ai^"'^^"+'^ = 2 - ^ = so that the Laplace exponent of 

the convolution semigroup z/^ defining the branching mechanism A^^^'^"^'''^^"'^ — A|j^"f ^'^equals 
— 1 = = /(I — e~^^)e~^ds. It is the semigroup of a compound Poisson process 
whose Levy measure is exponential. 

The immigration law (on IR+) is a Gamma distribution r(a,G^'^) = T{a, 1). It is the 
law of jCI^ and also of [vC^""^]" (it means the occupation field of the trace of Ca on Dn-i 
evaluated at n) for all n > 1. 

The conditional law of given jO^ is the convolution of the immigration law P(q;, 1) 
with 

Exercise 62 Alternatively, we can consider the integer valed process iV„(£i^"^) + l which 
is a Galton Watson process with immigration. In our example, we find the reproduction 
law Tr{n) — 2~'^~^for all n > (critical binary branching). 

Exercise 63 Show that more generally, if Cn,n+i = [iz-]"? for n > and ki — l,with 
Q < p < 1, we get all asymetric simple random walks. Show that A„ = ^fz^ andG^'^ — 1. 
Determine the distributions of the associated branching and Galton Watson process with 
immigration. 

If we consider the occupation field defined by the loops going through 1, we get a 
branching process without immigration: it is the classical relation between random walks 
local times and branching processes. 
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9 Loop erasure and spanning trees. 

Recall that an oriented link g is a pair of points {g^,g^) such that Cg = Cg~^g+ 7^ 0. 
Define -g = ig^,g~). 

Let fi^y be the measure induced by C on discrete self-avoiding paths between x and 

Another way to define a measure on discrete self avoiding paths from x io y from 
a measure on paths from x to y is loop erasure (see [13] ,[31], [H] and [2l])- In this 
context, the loops, which can be reduced to points, include holding times, and loop 
erasure produces a discrete path without holding times. 

We have the following: 

Proposition 64 The image of /i^'^ hy the loop erasure map 7 — * 7^'^ is fi^^ defined 
on self avoiding paths by HbUv) = f^""/ iv) dS§§F) = At^'^('7) det(G|{^}x{,,}) (Here {77} 
denotes the set of points in the path rj) and hy /i^|;(0) = SyGx,x 

Proof, li rj = (xi = x, X2, = ?/),and rjm = {x, ...,Xm), rn > 1 

00 

^^,y^^BE ^^^^ 5^[Pl^P:,/igf.(7^^ = ev) 

where denotes the bridge measure for the Markov chain killed as it hits x and 6 the 
natural shift on discrete paths. By recurrence, this clearly equals 



as 



for all m < n — 1. ■ 

Also, by Feynman-Kac formula, for any self-avoiding path 77: 

-<^'^>l|,..=,|/i^'nc^7) = ^^''l^^iy /iv) = det(G,)|{,|x{,}/x^/(r]) 

det(Gi ) 

^ det(Gx)lWxW (j^) 
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Therefore, recalling that by the results of section [573] conditionally on rj, Ci/Cl"^^^ and 
C\ are independent, we see that under /x"^'^, the conditional distribution of 7 given 
^BE = ^ ig ^j^g distribution of Ci — C\^^ i.e. the occupation field of the loops of Ci which 
intersect rj. 

More generally, it can be shown that 

Proposition 65 The conditional distribution of the network defined by the loops of 
7, given that ^y^^ = rj, is identical to the distribution of the network defined by Ci/C\^^ 
i.e. the loops of Ci which intersect rj. 

Proof. Recall the notation = det(G'). First an elementary calculation using ([9]) 
shows that /ig;^(e*'^T'^l{^sE=^}) equals 



Cx,X2Cxi,X-J,---Cxri-l,y U^V ^"'^ 

(Note the term e^^^ "^'^) can be replaced by ]!«( v")^"^^^)- 



Moreover, by the proof of the previous proposition, applied to the Markov chain de- 
fined by e' perturbed by we have also /i^;^(e'^^'^l|^BB=^|) = C^,^^ C'^^ • • -C^-i ,y ^ 1 Z 
Therefore, 



Moreover, by fl29l) and the properties of the Poisson processes. 

It follows that the joint distribution of the traversal numbers and the occupation field 
are identical for the set of erased loops and Ci/C\ . ■ 

The general study of loop erasure which is done in this chapter yields the following 
result when applied to a universal covering X . Let G be the Green function associated 
with the lift of the Markov chain. 
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Corollary 66 The image of jj,^'^ under the reduction map is given as follows: If c is a 
geodesic arc between x and y:fi^'^{{^,^^ = c}) = n^c^.c^+i det(G|{c}x{c}) 

Besides, if x and y are the endpoints of the lift ofc to a universal covering, ^^'^{{E^, C,^ = 
c}) = Gj.y 

Note this yields an interesting identity on the Green function G. 
Exercise 67 Check it in the special case treated in proposition l21\ 
Similarly one can define the image of by BE which is given by 

for T] = {xi, ...,Xn), and get the same results. 

Note that in particular, P^s(0) = V^{1 - Y.yPy) = f^xG'''''. 

Wilson's algorithm (see [23]) iterates this construction, starting with x's in arbitrary 
order. Each step of the algorithm reproduces the first step except it stops when it hits 
the already constructed tree of self avoiding paths. It provides a construction of a random 
spanning tree. Its law is a probability measure P|ji on the set STx,a of spanning trees 
of X rooted at the cemetery point A defined by the energy e. The weight attached to 
each oriented link g = {x,y) of X x X is the conductance and the weight attached to 
the link (x. A) is which we can also denote by C^^a- As the determinants simplify, the 
probability of a tree T is given by a simple formula: 

nTm = Zei[Ci (30) 

It is clearly independent of the ordering chosen initially. Now note that, since we get a 
probability 

^ ^ ^x,y '^x = 1 

Te5Tx,A {x,y)&T x,(x,A)eT 

or equivalently 

E n n ^--in^ 

Then, it follows that, for any e' for which conductances (including k') are positive only 
on links of e. 



11 pa: 11 pa; I TT V 7 



62 



and 

n §^ n (31) 

\{x,y)er ""^y x,{x,A)&r ^/ 

Note also that in the case of a graph (i.e. when all conductances are equal to 1), 
all spanning trees have the same probability. The expression of their cardinal as the 
determinant is Cayley's theorem (see for exemple p3|). 

The formula fIST]) shows a kind of duality between random spanning trees and Ci. It 
can be extended to Ck for any integer k if we consider the sum (in terms of number of 
transitions) of k independent spanning trees. 

Exercise 68 Show that more generally, for any tree T rooted in A, 

P|j.({T, T C T}) = det {G\{T}x{T}) Ylfer^i' {^} denoting the vertex set ofT. 



Corollary 69 The network defined by the random set of loops Cw constructed in this 
algorithm is independent of the random spanning tree, and independent of the ordering. 
It has the same distribution as the network defined by the loops of Ci. 

This result follows easily from proposition [651 

Remark 70 Note that proposition\6^ and its corollary can be made more precise with the 
help of remark [5M The splitting procedure used there with the help of an auxiliary inde- 
pendent set of Poisson Dirichlet variables allows to reconstruct the set of loops Ci/ 
by splitting the first erased loop in the proof of the proposition. Iterating the procedure 
we can successively reconstruct all sets c\^"^^ finally Ci/ . Then, by 

Wilson algorithm, we can reconstruct C\. 



Let us now consider the recurrent case. 

A probability is defined on the non oriented spanning trees by the conductances: 
P|y((T) is defined by the product of the conductances of the edges of T normalized by 
the sum of these products on all spanning trees. 

Note that any non oriented spanning tree of X along edges of E defines uniquely an 
oriented spanning tree Ia{T) if we choose a root A. The orientation is taken towards the 
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root which can be viewed as. a cemetery point. Then, if we consider the associated Markov 
chain killed as it hits A defined by the energy form e^'^^" , the previous construction yields 
a probability P|i^^ on spanning trees rooted at A which by fl30l) coincides with the image 
of Fgrp by /a- This implies in particular that the normalizing factor 2^.^^^^ is independent 
of the choice of A as it has to be equal to {^^^stx ^^{x yjeT C^x,y)~^- We denote it by Z^. 
This factor can also be expressed in terms of the recurrent Green operator G. Recall it is 
defined as a scalar product on measures of zero mass. The determinant of G is defined as 
the determinant of its matrix in any orthonormal basis of this hyperplane, with respect 
to the natural Euclidean scalar product. 

Recall that for any x ^ A, G{ex — £a) = — ' ^^j^-, + G^^^^'^e^- Therefore, for any 
y^A, {By - G{e, - e^)) = [G^'^^-'y. 

The determinant of the matrix [G^^J'''], equal to Z^^ is therefore also the determinant 
of G in the basis {5x — (5a, x 7^ A} which is not orthonormal with repect to the natural eu- 
chdean scalar product. An easy calculation shows it equals det((5y — 5a, — '5a)]r|xi ,x,y ^ 
A) det(G) = |X| det(G). Note also that if we set ax^iT) = n(x,j/)G/:.„(r) ^t<^sTx '^^oi'^) 
is proportional to X^^ as Xq varies in X. 

Exercise 71 Prove that if we set a^^iT) = n(x,3/)e/,,)(r) -^y ' Y^reSTx '^^oi'^) propor- 
tional to X^^ as Xq varies in X . More precisely, it equals KXxq, with K = — This 
fact is known as the matrix- tree theorem (JEE/)- 

Exercise 72 Check directly, that Zg^^^^^ is independent of the choice of xq. 

Let us come back briefiy to the transient case by choosing some root xq. As by the 
strong Markov property, = fyiT^ < oo)V^, we have = M = ^yiT^ < 00), and 
therefore 

nAix,y) G T) = P.(7r = y) = V:Py^Fy{Tx = 00) = Gx,yG^'^il - ^). 
Directly from the above, we recover Kirchhoff's theorem: 



P|^(±(x, y)eT) = Gx,y[G^'^{l - ^) + G^'^l - ^)] 

= Cx,y{G'''^ + G^'y - 2G'^'y) = Gx,yK'''y^''''''y^ 
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with the notations of section 12. 5[ and this is clearly independent of the choice of the 
root. 

Exercise 73 Give an alternative proof of Kirchhoff's theorem by using l[3l\) . taking C'^y = 
sC^^y and Q„ = Cu,v for {u,v} ^ {x,y}. 

More generally we can show the transfer current theorem (see for exemple [22], [23]): 

Theorem 74 P|t(±^i, ... ± G T) = (IlJCgJ det(ir«-«^ 1 < i,j < k) 

Note this determinant does not depend on the orientation of the links. 

Proof. Note first that if T is a spanning tree rooted in xq = A and = Xj), 1 < 

i < 1^1 ~ 1 cire its oriented edges, the measures — form another basis of the 
euclidean hyperplane of signed measures with zero charge, which has the same determinant 
as the basis — 5xo- 

Therefore, is also the determinant of the matrix of G in this basis, i.e. 

Zl = det{K^^'^^ 1 < J < 1^1 - 1) 

and F'st{T) = (]ll^'"^ det(i^«-«^ l<i,J < \X\ - 1) 
= det(^/C^ir«-«^Y4^ 1 <hJ < \X\ - 1). 

Recall that ^JU^^K^^'^^ ^JC^. = |n|Q;^^, ^ , where 11 denotes the projection on the 
space of differentials and that aT-!^ = ±— ^ if {x,y) = ±{ri) and = elsewhere. 

To finish the proof of the theorem, it is helpful to use the exterior algebra. Note 
first that for any ONB ei, ...e|x|-i of the space of differentials, Ila^ = Yl 

P|t(T) = det((a^Je,»2 = (^al A ... A a^^^^Je, A ... A e|x|-i)^|_,|_^^ • Therefore 

P|^(ei, , e T) = E,,,„...„,,_, («6 A ...al A a;^^ A Je^ A ... A e|x|.i)^,, 

where the sum is on all edges rj^+i, ...?7|x|-i completing ^i, , ...,^k into a spanning tree. It 
can be extended to all systems of |X| — 1 — edges as all the additional term vanish. 
Indeed, the projection on an exterior product of differentials of an exterior product of 
vanishes as soon as they form a loop. Hence the expression above equals: 

En<i....<i. R A ...a|Je,, A ... A e,,)' ^_ = det(^ir«''«^V^ 1 < z, j < k) m 
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It follows that given any function g on non oriented links, 
E|^(e-E,.r.(0) = E|^(n(l + (e-^^«^ - l)her) 

= E E n(^"'^^'^ - i)P|t(±6, ... ± a e T) 

= E E n^^""^^"^ ~ ^) det(ir«-^^ 1 < ^, J < fc) 

= 5^ Tr((Mc(e-.-i)ir)^') = det(J + KMcie~.-i)) 
and we have proved the following 

Proposition 75 E|^(e-See-^(«)) = det(/ - M ^y^^^-^^K M ^^^^^—^) 

Here determinants are taken on matrices indexed by E. 
This is an exemple of the Fermi point processes discussed in [35]. 

As a consequence, for any spanning tree T, if ttt denotes Mi^^^ (the multiplication by 
the indicator function of T), it follows from the above, by letting g be mli^c}, m ^ oo 
that P|t(T) = det((/ - KMc){I - vtt) + vtt) = det((/ - KMc)t-xt-) 

c 

Note also that if e' is another energy form on the same graph, IE|T(n(x j/)eT cf^) ~ 
det{I-M^z^KM 



On the other hand, from (I3T1) . it also equals (X] 
finally 



20 

^ SO that 



ll = det(/ - M 



c-c 



V c V c 

Note that indicators of distinct individual edges are negatively correlated. More gen- 
erally: 

Theorem 76 (Negative association) Given any sets disjoint of edges Ei and E2, Fgrp[Ei[J 
C T) < P|r(^i C T)F'st{E2 C T) 



Proof. Denote by K'^{i,j) the restriction of K"^ 
E, X E2. Then, nr^E.^E.^r) _ det(x#i 

Ir c 



K*{l,l)-hK*{l,2)K*{2,2)~h . 



Aet{K*{2,2)) Aet{K*(2,2)) 



/Qir^'^V^, i.V ^ E) to 
I F 
F* I 



det( 



with F 
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Note that the matrix F is of the form (cj, fj) where Cj and fj are systems of ortho nor- 



mal vectors in 
determinant of 
M^i and 



l^'^ and for the scalar product defined by fC*. It follows that the 
I F 

is the square of the sine of the angle between the two subspaces 



10 Reflection positivity 

In this section, we assume there exists a partition of X: X = X+ U X^, X+ n X^ = 
and an involution p on X such that: 

a) e is p-invariant. 

b) p exchanges and X~. 

c) The X+ X X+ matrix C^^y = Cx,p(y), is nonnegative definite. 
Then the following holds: 

Theorem 77 ij For any positive integer d and square integrable function $ in a{Cd , x G 

x+) vor(xi3,x,?/eX+), 

E($(£,)$(p(£,)) > 
a) For any square integrable function S of the free field (j) restricted to X^ , 

E^(E((/.)E(p(0)) > 

Hi) For any set of edges {^i} in X^ x X+ the matrix, 

K.,j = P5t(6 e T, Pi, G T) - ¥sT{^^ e T)FsT{ij e T) 

is nonpositive definite 

Proof. The property ii) is well known in a slightly different context: Cf for exemple 



|lUj ] and their references. Reflection positivity is a keystone in the bridge between 
statistical and quantum mechanics. 

To prove i) , we use the fact that the cx-algebra is generated by the algebra of random 
variables of the form $ = ^ ^^■5^'^''^^ with = C and Uj = on X" x X", C(<^i) < C, 
Afe) = A on X- and A^-^^) > A. 
Then 
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y with A(^^.«) 



c 



We have to prove this is non negative. It is enough to prove it for d = 1, as the 
Hadamard product of two nonnegative definite Hermitian matrices is nonnegative definite. 

Let us first assume that the nonnegative definite matrix is positivedefinite. We 
will see the general case can be reduced to this one. 

Now note that 2g.^p(^eq)-e,u).j-p{u)q) is the inverse of the determinant of a positive definite 



matrix of the form: D(j, q) — 



and C^'''^ = D^^W^"^ with - , 

It is enough to show that det(D(j, can be expanded in series of products ^ qn{j)Qn{k) 
with \Qnij)f < oo. 



A{j) -C±'^'^' 
_C±:hj A{q)* 
1 



with = A^^(5, 



u,v 



u,p(u.) 



As 



D{3. q) 









-I 



the — a-power of this determinant can be written 

det(L'(^))-Met(:D^'V'det(F(j)) det(F(?)*) det(7 ■ 
with F{j) = D(^^A{j)-^D(^\ or more simply: 





F{qy 



m 





F(j) = det(A(j))-Met(A(?)*)-Met(7 



F{j) 
F{qy 

Note that A{j)~^ is also the Green function of the restriction to X+ of the Markov 
chain associated with ej, twisted by Uj. Therefore A(j)~^D^^^ — ■§^[C^''^'^]''^ F(j)[C^]^ 
is the balayage kernel on X~ defined by this Markov chain with an additional phase under 
the expectation produced by uj. It is therefore clear from Frobenius theorem that the 
eigenvalues of the matrices A{j)'^D^^^ and F{j) = ^^A{j)~^D''^'' are of modulus less than 



one and it follows that 



F{j) _ I F{q)* 
F{qy J ~ [ / J [ F{j) 

If X+ has only one point, (1 - F{j)F{q)*)-^ = J^FU^Hq) 
conclude. Let us now treat the general case. 



is a contraction, 
which allows to 
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For any (n,m) matrix A^, and k = {ki,...km) e N"", / = {h,...ln) e N", let A^^'^''^ 
denote the {\k\ , matrix obtained from by repeting ki times each line i; then Ij times 
each column j. 

We use the expansion 



det(/ - M) 



-1 



valid for any strict contraction M (Cf |10] and [S]). 

Note that if X has 2d points, if we denote (fci, ...k2d) by A;^), with k^ = {ki, ...kd) 



and k 



.k 



2d 



F{j) 
F{q)* 



{fc,A;} 







[F(^qy]{k-,k+} 

But the all terms in the permanent of a (2n, 2?t,) matrix of the form 









vanish 



A 

B* 

unless the submatrices A and i? are square matrices (not necessarily of equal ranks). 
Hence in our case, we necessary have = \k~\, so that, A and B are (n, n) matrices. 

Then, the non zero terms in the permanent come from permutations exchanging 
{l,2..n} and {n + l,...2n}, which therefore decompose into a pair of permutations of 
{l,2..n}. Therefore : 



Per{ 



A 

B* 



Per{A)Per{B*) 



which concludes the proof in the positive definite 

Per (5*) = Eres. Hi = Er,s. Hi ^^W. = P^)- 

To treat the general nonnegative case, we can use use a passage to the hmit or alter- 
natively, the proposition [60] (or more precisely its extension including a current) to reduce 
the sets and X~ to the support of C^. 

To prove iii) let us first show the assumptions imply that the X+ x X+ matrix G^, 

r A — 

Gx,p{y) is also nonnegative definite. Let us write G in the form . 
A = Mx-C. Then 



with 



G 









A~h 



-A~iC^A-' 



-A-W^A~2 
I 



-, -1 


' " 




_ A-^ _ 



A-2C^A- '2 is non negative definite and as before, we can check it is a contraction 
since A~^C^ is a balayage kernel. 
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Note that if a symmetric nonnegative definite matrix K has eigenvalues /Xj, the eigen- 

^ D E 



values of the symmetric matrix E defined by 



I -K 
-K I 



E D 



are easily seen 



(Exercise) to be -yz^- Taking K = A ^C^A 2^ it follows that the symmetric matrix E 
, (and in our particular case = A-2EA-2) is nonnegative definite. 
To finish the proof, let us take S of the form ^ A^e^'^'^^^ Then 

E^(S(0)E(p(0)) = ^A,-A,E^(e<*'>^^>+<'^'^('^')>) = ^ A,-eHx..G++x.)^^eHx..G++x,)eM>^^'^*x.) 

and as is positive definite, we can conclude since e^^^^'*^ ^'') = Eu,(e^"''-^^^e^"''-^''^^), 
w denoting the Gaussian field on with covariance G^. 

To prove iii), note that Ki^j = —{Gf^Y- Then, using again the Gaussian vector w, 
and the Wick squares of its components; {Gf.jy = E(: wf :: w'j :). ■ 

Remark 78 a) If Uj are unitary representations such that Uj'^ is the identity outside 
X X^, i) can he extended to variables of the form ^ XjB^^^-^'^^ and to the a-field they 
generate. 

b)In the case where a is a half integer, by remarklj^l the reflection positivity of the 

X 

free field ii), implies i) holds also for any half integer a provided that $ G a{Ca ,x G 
X+) Va(ivS + iV("j,x,i/GX+). 

Exercise 79 Prove the above remark. 

Remark 80 // there exists a partition of X : X = X~^ U X~ U and an involution p 
on X such that: 

a) e and X^ are p-invariant. 

b) p(X±) = X^ 

c) X+ and X" are disconnected. 

Then the assumptions of the previous theorem are satisfied for the trace on X^ U X^ . 

Moreover, if X^ x X^ does not contain any edge of the graph, the assertion i) of 
theorem 1771 hold for the non disjoint sets X^ U and X~ U X". More precisely, i), 
holds for $ in 0" ,x,y e It is enough to apply the theorem 

to the graph obtained by duplication of each point xq in X^ into {xq,Xq) and let all 
conductances ^- equal to a constant increasing to infinity. 

Physical Hilbert space and time shift: We will now work under the assumptions 
of the previous remark, namely, without assuming that X = X~^ U X~ . 
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Let TC^ be the space of square integrable functions in cr((>Ci , x e X+) V a{Ni^y,x, y e 
X+), equipped with the Hermitian scalar product = E($(>Ci)^'(p(>Ci)) Note that 

($,$)^<E($$(A))- 

Let A/" be the subspace G 7i:+, E(^(£i)\l>(p(£i)) = 0} and H the closure (for the 

topology induced by this scalar product) of the quotient space Ti^ /M (which can be 
called the physical Hilbert space). We denote $~ the equivalence class of $. Note that 
$ — is equipped with the hermitian scalar product defined unambiguously by 

Assume X is of the form x Z (space xtime) and let 9 be the natural time shift. 
We assume 9 preserves e, i.e. that conductances and k are ^-invariant. We define p by 
p{xo,n) = {xo, —n) and assume e is p-invariant. Note that 9{X~^) C and p9 — 9~^p. 
The transformations p and 9 induce a transformations on loops that preserves p, and 9 
induces a linear transformation of Ti"*". Moreover, given any F in A/", F o 9 & J\f, as 

{Fo9,Fo 9)^ is nonegative and equals E{F o 9{£i)F o ^(p(>Ci)) = E{F{9{£i))F{p o 
ri(A)) = E(F(g^(£i) )F(p(£0) 

= (F o F)„ < y/{F o9\F o9^)^{F,F)^ which vanishes. 

Proposition 81 T/?,(;re exist a self adjoint contraction of Ti, we will denote U.^^^ such 
that [$0^]- = n(^)($~) 

Proof. The existence of n^^-' follows from the last obsevation made above. As 9 pre- 
serves Fc^, it follows from the identity p9 = 9-^p that {F o 9, G)^ = E{F{9{Ci))G{p{Ci) = 
E{F{Ci)G{p o 9-\Ci)) = E{F{Ci)G{9 o p{C^)) ^{Fo9, G)^. Therefore, is self ad- 
joint on T-C^ /M. To prove that it is a contraction, it is enough to show that {F o9,F o 9)-^ < 
{F,F)^ for all F en+. 

But as shown above, {Fo9,Fo 9)^ = (F o 9^, F)^ < s/iF o 9^, F o 9^)^ {F,F)^ 

By recursion, it follows that: 

{Fo9,Fo 9)^ = {F o 9\ F)^ <{Fo9'\Fo 9'%" {F, F)]^''" 

As {F o 9^'\F o 9^"^^" < (E(F(£i)F(£i))2"" converges to 1 as n ^ oo, the inequal- 
ity follows. ■ 

For all n G Z, the symmetry p*^") = 9~'^p9^ allows to define spaces isometric to 
H. These isometrics can be denoted by the shift For n > m, jn,m = ^™[n(^)]"-"'^-" 
is a contraction from TY^") into Tt^'^^ . 
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11 The case of general Markov processes 



We now explain briefly how some of the above results will be extended to a symmetric 
Markov process on an infinite space X. The construction of the loop measure as well 
as a lot of computations can be performed quite generally, using Markov processes or 
Dirichlet space theory (Cf for example ^). It works as soon as the bridge or excursion 
measures can be properly defined. The semigroup should have a locally integrable 
kernel y). 

Let us consider more closely the occupation field /. The extension is rather straight- 
forward when points are not polar. We can start with a Dirichlet space of continuous 
functions and a measure m such that there is a mass gap. Let Pt denote the associated 
Feller semigroup. Then the Green function is well defined as the mutual energy of the 
Dirac measures 5x and 5y which have finite energy. It is the covariance function of a 
Gaussian free field which will be associated to the field C\ of local times of the 

2 

Poisson process of random loops whose intensity is given by the loop measure defined 
by the semigroup Pt- This will apply to examples related to one- dimensional Brownian 
motion or to Markov chains on countable spaces. 

When we consider Brownian motion on the half line, the associated occupation field 
Ca is a continuous branching process with immigration, as in the simple random walk 
case. 

When points are polar, one needs to be more careful. We will consider only the case 
of the two and three dimensional Brownian motion in a bounded domain D killed at the 
boundary, i.e. associated with the classical energy with Dirichlet boundary condition. 
The Green function does not induce a trace class operator but it is still Hilbert- Schmidt 
which allows us to define renormalized determinants det2 (Cf |33j). 

If A is a symmetric Hilbert Schmidt operator, det2(/ + A) is defined as ^(1 + Aj)e~'^' 
where Aj are the eigenvalues of A. 

The Gaussian field (called free field) whose covariance function is the Green function 
is now a generalized field: Generahzed fields are not defined pointwise but have to be 
smeared by a test function /. Still 0(/) is often denoted / (l){x)f{x)dx. 

Wick powers : cf)^ : oi the free field can be defined as a generalized field by approxi- 
mation as soon as the 2n-th power of the Green function, G{x,yY^ is locally integrable 
(Cf [31]). This is the case for all n for Brownian motion in dimension two, as the Green 
function has only a logarithmic singularity on the diagonal, and for n = 2 in dimension 
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three as the singularity is of the order of -pc^- More precisely, taking for example 7T^{dy) 
to be the normalized area measure on the sphere of radius e around x, 0(vr^) is a Gaus- 
sian field with variance = / G{z, z')TT'^{dz)7!'^{dz'). Its Wick powers are defined with 
Hermite polynomials as we did previously: 

: 0(vr|)" := (cr^)tif„(^^^). Then one can see that, / f{x) : 0(vr|)" : dx converges in 



for any bounded continuous function / with compact support towards a limit called 
the n-th Wick power of the free field evaluated on / and denoted : 0" : (/). Moreover, 



In these cases, we can extend the statement of theorem [36] to the renormalized occu- 
pation field £| and the Wick square : 0^ : of the free field. 



Let us explain this in more detail in the Brownian motion case. Let D be an open 
subset of M'^ such that the Brownian motion killed at the boundary of D is transient and 
has a Green function. Let pt{x,y) be its transition density and G{x,y) = pt{x,y)dt 
the associated Green function. The loop measure /i was defined in [I5j as 



where P^'^ denotes the (non normalized) excursion measure of duration t such that if 
< ti < ...th < t, 

^t'^'Uih) e dxi,...,^{th) e dxh) = pt^{x,Xi)pt2^t^ixi,X2) pt-th{xh,x)dxi...dxh 

(the mass of P^'"^ is Pt{x, x)). Note that is a priori defined on based loops but it is easily 
seen to be shift-invariant. 

For any loop / indexed by [0 T(/)], define the measure / = Si(^s)ds: for any Borel 



{x,y)f{x)h{y)dxdy. 




set A, T{A) = J^^'^ UiQds. 



Lemma 82 For any non-negative function f 



Ai((^/) ) = ('^-l)' / G{Xi,X2)fix2)G{x2,X3)f{x3)...G{Xn,Xi)f{Xi)YldXi 
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Proof. Prom the definition of /i and I, //(^/, /y ) equals: 
i! / / ^f{xi)...f{xn)ptAx,Xi)....pt-t„{xn,x)Y\dtidxidtdx 

J J {0<ti...<tn<t} 




^f{xi)...f{Xn)pt2-ti{xx,X2)...Pu+t-tSXn,Xi)\YdUdXidt 

{0<tl...<tn<t} ^ 

Performing the change of variables V2 — t2 — ti, ..Vn — tn — tn-i,vi — ti + t — tn, and 
V — ti, we obtain: 

n\ I f{Xi)...f{Xn)Pv2{Xl,X2)....Pvi{Xn,Xi)Y\dVidXidv 

J{0<v<vi ,0<Vi} Vi + ...+Vn 

f{xi)...f{Xn)Pv2{xi,X2) Pv,{Xn,Xi)Y\ dVidXi 

J{0<Vi}'Vl + ■■■ + Vn 

= (n-l)!/ /(.,).../Wp„(.....)...p„(.....)n'"'*. 

J{0<v,} 

(as we get the same formula with any Vi instead of Vi) 

/n 
G{xi, X2)f{x2)G{x2, X3)f{x3)...G{Xn, Xi)f{xi) dXi- 
1 

■ 

One can define in a similar way the analogous of multiple local times, and get for their 
integrals with respect to a formula analogous to the one obtained in the discrete case. 

Let G denote the operator on LF'{D,dx) defined by G. Let / be a non- negative 
continuous function with compact support in D. 

Note that (l^ is //-integrable only in dimension one as then, G is locally trace class. 

In that case, using for all x an approximation of the Dirac measure at x, local times P 
can be defined in such a way that <^/, = / l^f{x)dx. 

(^l, is yU-squarc intcgrable in dimensions one, two and three, as G is Hilbert-Schmidt 
if D is bounded, since J Jj^^jyG{x,y)^dxdy < oo, and otherwise locally Hilbert-Schmidt. 

N.B.: Considering distributions x such that / J {G{x,y)'^x{dx)x{dy) < oo, we could 
see that (l,x) can be defined by approximation as a square integrable variable and 
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l^i(tx) ) = J{G{x,yfx{dx)x{dy). 



Let 2; be a complex number such that Re(2;) > 0. 

Note also that e + z (l, f^j — 1 is bounded by ^ and expands as an 

ahernating series f (" {l /))"> with |e-^(^^) - 1 - Ef (- {l /))"| < 
Then, for \z\ small enough., it follows from the above lemma that 

\ I ^ n 

As M^GM^ is Hilbert-Schmidt det2(/ + zM^GM^) is well defined and the second 
member writes -log(det2(/ + zM^GM^)). 
Then the identity 

^{e-<^'S) + JT,f)-l)^- log(det 2(/ + zM^GM^)). 



extends, as both sides arc analytic as locally uniform limits of analytic functions, to all 
complex values with positive real part. 

The renormalized occupation field is defined as the compensated sum of all I 
in Ca (formally, Ca — Ca — j j^^^^ 5i^ds^{dl)) By a standard argument used for the 
construction of Levy processes. 



^a, fj = Jim (^Ca,e, f 

with by definition 

The convergence holds a.s. and in L^. as 

E(( E,e£. (1{T>.} r !{ls)ds)-a^{\T>e} Jo f{ls)ds)r) = a /(1{t>.} f{is)dsfl^{dl) 
and E(^>Cq,, ) = Tr{{M^GM^Y). Note that if we fix /, a can be considered as a 



time parameter and (^Ca,e, fj as Levy processes with discrete positive jumps approximat- 
ing a Levy process with positive jumps (^Ca,f^- The Levy exponent /i(l{T>£}(e~('''^) + 
l,f) — 1)) of (Ca,£,f/) converges towards the Levy exponent of (Ca,f/) which is 
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^((e-('^/) ^(Tj'j- 1)) and, from the identity E(e-(^"'^)) = e'^'^^' + get 
the 

Theorem 83 Assume d < 3. Denoting the compensated sum of all I in Ca, we have 
E(e'(^'^)) = det2(/ + M^GM^))-^ 

Moreover e"^^"-''-^) converges a.s. and in towards e~^'~'°"-^\ 

Considering distributions of finite energy X (i.e. such that J{G{x,y)^x{d^)x{dy) < oo), 
we can see that (^Ca, xj can be defined by approximation as Yimx-^ooii^/^a^ ^G\x^) and 

E((Z;,x)') =a j {G{x,y))\{dx)x{dy). 
Speciahzing to a = |, /c being any positive integer we have: 

Corollary 84 The renormalized occupation field Ck and the Wick square | : Yli 'Pi '■ 
have the same distribution. 



If G is a conformal map from D onto O(-D), it follows from the conformal invariance of 
the Brownian trajectories that a similar property holds for the Brownian" loop soup"(Cf 
[15]). More precisely, if c(x) = J acohiarixiQ) and, given a loop /, if T'^il) denotes the 
reparametrized loop Ir^, with j^" c{lu)du = s, the configuration QT'^{Ca) is a Brownian 
loop soup of intensity parameter a on O(-D). Then we have the following: 

Proposition 85 Q{cCa) is the renormalized occupation field on Q{D). 

Proof. We have to show that the compensated sum is the same if we perform it after 

or before the time change. For this it is enough to check that 

H[^i'^{TT>v}'^{T<e} / fi%)ds-a il{rT>v}kT<e} / f i7s)ds) ^{d^f)]'^) 
= a j {l{TT>v}'^{T<e} f{ls)dsffi{dj) 
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and 

H[^i'^{T>e}UT<v / f{ls)ds-a j {l{T>e}^TT<n / 



f{js)ds)fi{dj)] 



2\ 



T 

2, 



a j {l{T>e}lTT<v j fils)ds) IJ,{d-f) 

converge to zero as e and t] go to zero. It follows from the fact that: 

j[kT<e} j^^ f{ls)dsfKdl) 

and 



converge to 0. The second follows easily from the first if c is bounded away from zero. 
We can always consider the "loop soups" in an increasing sequence of relatively compact 
open subsets of D to reduce the general case to that situation. ■ 

As in the discrete case (see corollary [32]) , we can compute product expectations. In 
dimensions one and two, for fj continuous functions with compact support in D: 

E( (ZaJi)... (Ca, fk)) = j Per° 1 < /, m < fc) J] fj{xj)dx,. (32) 

12 Renormalized powers 



In dimension one, as in the discrete case, powers of the occupation field can be viewed 
as integrated self intersection local times. In dimension two, renormalized powers of the 
occupation field, also called renormalized self intersections local times can be defined as 
follows: 

Theorem 86 Assume d = 2. Let 7c^{dy) be the normalized arclength on the circle of 
radius e around x, and set = J G{y, z)n^{dy)TT'^{dz) . Then, J f{x)Q^''^^{(^Ca,T^f^)dx 
converges in for any hounded continuous function f with compact support towards a 
limit denoted (^C^, and 

E((£^, /) (ZI, h)) = 5,,,£^(f^±l):^|fi±^ / G^'ix, y)f{x)h{y)dxdy. 
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Proof. The idea of the proof can be understood by trying to prove that 

^iij /(^)Qfc''^'" ((^-^a' z^)*^^)^) remains bounded as e decreses to zero. One will ex- 
pand this expression in terms of sums of integrals of product of Green functions and check 
that the combinatorial identities fl2T|) imply the cancelation of the logarithmic divergences. 

This is done by showing (as done below in the proof of the theorem) that one can mod- 
ify slightly the products of Green functions appearing in K{Q'^''^^ {(^'C.a, '^e'^)Qk'"^ {(^'^a, ^rf ^)) 
to replace them by products of the form G{x,yy{a^y{a^)^ . The cancelation of terms 
containing and/or then follows directly from the combinatorial indentities. 

Let us now prove the theorem. Consider first, for any xiX2---Xn, s small enough and 
e < El, ...En < ^E, with Ei = Ej if Xi = Xj, an expression of the form: 



A 



— 1 7~*^i 



in which we define rrii as sup(/i,a;i+fe = Xi). 

In the integral term, we first replace progressively G{yi-i,yi) by G{xi-.i,Xi) whenever 
7^ Xi, using triangle, then Schwartz inequalities, to get an upper bound of the absolute 
value of the difference made by this substitution in terms of a sum A' of expressions of 
the form 

l[G{xi, ^ j{G{y,,y2)-G{x,,X2))H'!l{dyMe'.{dy2) j \{G^{yk,yk+i)\{7^et{dyk). 
The expression obtained after these substitutions can be written 

W = JJ G{xi-i,x,) j G{yi,y2)...G{ym,_^,ymMl{dyi)-<:{dymy) 

and we see the integral terms could be replaced by (^cr^i'^'^i if Q was translation invariant. 
But as the distance between x and y tends to 0, G{x,y) is equivalent to Go{x,y) = 
^\og{\\x — y\\) and moreover, G{x,y) = Go{x,y) — H^\x,dz)Go{z,y), H^" denoting the 
Poisson kernel on the boundary of Z^. As our points lie in a compact inside D, it follows 
that for some constant C, for — x|| < |/(G'(?/i, 7/2)71^(^1/2) — erf | < Ge. 
Hence, the difference A" between W and Hixi-i/xi '^^^ be bounded by 

eW , where W is an expression similar to W . 
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To get a good upper bound on A, using the previous observations, by repeated applica- 
tions of Holder inequality, it is enough to show that for e small enough and £ < £i, £2 < 2s, 
(with C and C denoting various constants): 

1) /(Gd/i.ya) -G(xi,,X2)2<;(rfyi)<|(rfi/2) 

2) / G{yi,y2)''Tr^{dyi)n^{dy2) < C |log(£)|'' and more generally 

3) /G'(?/i,y2)'</(ci?/i)<^^(dy2) <C|log(£)|'= 

As the main contributions come from the singularities of G, they follow from the 
following simple inequalities: 

1') 

j |log(£2 + 2Re cos{9) + R^) - \og{R)f d9 
= j \log{{e/Rf + 2{e/R) cos{9) + 1) |' d9 < C{{e^n>v^y + \og'{R/e)l^j,^^,yy) 

(considering separately the cases where ^ is large or small) 

2') / |log(£2(2 + 2 cos(^))) l'' de<C \\og{e) f 

3') / |log((£i cos(^i) + £2 cos(^2) + rf + {si sm{9,) + £2 sin(^2))'|' d9id92 < C(|log(£)|)^ 
It can be proved by observing that for r < ei + 62, we have near the singularities 
(i.e. the values 6i{r) and 62{r) for which the expression under the log vanishes) 
to evaluate integrals bounded by C jQ{—log{eu))''du < C"(— log(£))'^ for e small 
enough. 

Let us now show that for £ < £1, £2 < 2£, we have, for some integer N'^''^ 



< Clog(£)^'.'=(v^ + G'(x,y)^^l{||,_,||<^). (33) 
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Indeed, developing the polynomials and using formula (15^ we can express this expec- 
tation as a linear combination of integrals under JJ^ Tr^^^lfia^j) ]^^l2('^%) products of 

i j 

G{xi, Hi/), G{xi, Xj) and G{yj, yji) as we did in the discrete case. If we replace each (^(xj, Uj) 
by G{x, y), each G{xi, Xi') by cr^^ and each G{yj, yj') by a^^, we can use the combinatorial 

identity (I2T]) to get the value Si kG(x, yY^ — . Then, the above results 

k\ 

allow us to bound the error made by this replacement. 

The bound ( l33i) is uniform in (x, only away from the diagonal as G{x,y) can be 
arbitrarily large but we conclude from it that for any bounded integrable / and h, 

J (E(g^'-' {(Z, 7r:,))Qr" {{Z, T^l))) - kkGix, yy^^:::^^±^^)f(x)hiy)dxdy 

(as / / ^(x, can be bounded by Ge^, for example). 

Taking e„ = 2~", it is then straightforward to check that / f{x)Q'k"'' i.(^'^a,''^s„'^)dx 
is a Cauchy sequence in L^. The theorem follows. ■ 

Specializing to a = |, fc being any positive integer as before. Wick powers of Yl^=i 0^ 
are associated with self intersection local times of the loops. More precisely, we have: 

Proposition 87 The renormalized self intersection local times Ol and the Wick powers 

2 

2^ ■ (Xli • have the same joint distribution. 

The proof is similar to the one given in [19] and also to the proof of the above theorem, 
but simpler. It is just a calculation of the L^-norm of 

(0^r :(x)-gH(: 0^«))]/(x)dx 

which converges to zero with e. 

Final remarks: 

a) These generalized fields have two fundamental properties: 
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Firstly they are local fields (or more precisely local functionals of the field in the 
sense that their values on functions supported in an open set D depend only on the 
trace of the loops on D. 

Secondly, noting we could use different regularizations to define the action of 
a conformal transformation B on these fields is given by the k-th power of the 
conformal factor c = Jacobian{Q). More precisely, G(c'^£^) is the renormalized 
k-th power of the occupation field in Q{D). 

b) It should be possible to derive from the above remark the existence of exponential 
moments and introduce non trivial local interactions as in the constructive field 
theory derived from the free field (Cf |34j). 

c) Let us also briefly consider currents. We will restrict our attention to the one and two 
dimensional Brownian case, X being an open subset of the line or plane. Currents 
can be defined by vector fields, with compact support. 

Then, if we now denote by the complex valued free field (its real and imaginary 
parts being two independent copies of the free field), u and f^{(f)d^(f)—(f)d^(f))dx are 
well defined square integrable variables in dimension 1 (it can be checked easily by 
Fourier series). The distribution of the centered occupation field of the loop process 
"twisted" by the complex exponential exp(^^g^^ J^ia; + |/(||a;||^)) appears to be the 
same as the distribution of the field : 00 : "twisted" by the complex exponential 
exp(/_^(0c'^0-0c'^0)rfa;) (Cf[20]). 

In dimension 2, logarithmic divergences occur. 

d) There is a lot of related investigations. The extension of the properties proved here 
in the finite framework has still to be completed, though the relation with spanning 
trees should follow from the remarkable results obtained on SLE processes, especially 
[17] . Note finally that other essential relations between SLE processes, loops and 
free fields appear in [12], [32] and [6j. 
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